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Abstract. We study integrals of completely integrable quantum systems associated with 
classical root systems. We review integrals of the systems invariant under the corresponding 
Weyl group and as their limits we construct enough integrals of the non-invariant systems, 
which include systems whose complete integrability will be first established in this paper. 
We also present a conjecture claiming that the quantum systems with enough integrals given 
in this note coincide with the systems that have the integrals with constant principal symbols 
corresponding to the homogeneous generators of the i?„-invariants. We review conditions 
supporting the conjecture and give a new condition assuring it. 
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1 Introduction 

A Schrodinger operator 



1 n 

=1 3 



1.1 



with the potential function R(x) of n variables x = (x\, . . . , x n ) is called completely integrable if 
there exist n differential operators P±, . . . ,P n such that 



[Pi,Pj] = (l<i<j<n), 
P€C[P l ,...,P n ], 

Pi, . . . ,P n are algebraically independent. 



(1.2) 



In this paper, we explicitly construct the integrals P\,...,P n for completely integrable potential 
functions R(x) of the form 



R(x) = ^2 ( u ij( x i- x j)+ u tj( X i + X j)) + ^2 V k( x k) 
l<«<j<" fc=l 



(1.3) 



*This paper is a contribution to the Vadim Kuznetsov Memorial Issue 'Integrable Systems and Related Topics'. 
The full collection is available at |http://www.emis.de/journals/SIGMA/kuznetsov.html| 
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appearing in other papers. The Schrodinger operators with these commuting differential oper- 
ators treated in this paper include Calogero-Moser-Sutherland systems (cf. [5] |2"2" 1 |2"71 I2"5| 155]). 
Heckman-Opdam's hypergeometric systems (cf. |34j for type A n -\, [11] in general), their exten- 
sions (cf. [H [HJ HH [23j 1231 [251 EQ] ) and finite Toda lattices corresponding to (extended) Dynkin 
diagrams for classical root systems (cf. [21 [9j [10], [T71 [261 E31 EH]) arid those with boundary 
conditions (cf. PElttHKIlEOlEIlilEO]). 

Put dj = d/dxj for simplicity. We denote by cr(Q) the principal symbol of a differential 
operator of Q. For example, a(P) = — (l/2)(£f H 

We note that [3Qj proves that the potential function is of the form (jl.3p if 

*(flfc) = E for k = l,...,n. (1.4) 

l<h<— <jk<n 

In this case we say that R(x) is an integrable potential function of type B n or of the classical type. 
Moreover when R{x) is symmetric with respect to the coordinate (x\, . . . , x n ) and invariant under 
the coordinate transformation (x±, X2, ■ ■ ■ , x n ) t— > (— xi, a?2, . . . , x n ), then i?(x) is determined 
by [32] for n > 3 and by [23] for n = 2 and are calculated by [29j . 

Classifications of the integrable potential functions under certain conditions are given in 
[231 EH EU EOj [37l [40] etc. In Section [9j we review them and we present Conjecture which claims 
that the potential functions given in this note exhaust those of the completely integrable systems 
satisfying (jl.4h . We also give a new condition which assures Conjecture. 

If V/s = for k = 1, . . . , n, we can expect cr(Pk) = Yl ^ ' ' ' f° r & = 1; • ■ ■ > n ~ 1 

i<n<---<jfc<™ 

and cr(Pn) = i\£,2 ■ ■ -£,n and we say the integrable potential function is of type D n . If Vk = 
and = for k = 1, . . . , n and 1 < i < j < n, we can expect Pi = d\ + ■ ■ ■ + d n , c(-Pfc) = 
Cj 1 • • ■ <v/ fc for k = 2, . . . , n and we say that the integrable potential function is of 

l<jl<-<jk<n 

type A n _\. Note that the integrable potential function of type A„_i or D n is of type B n . 
The elliptic potential function of type A n _\ with 

= Cp(t; 2wi, 2wz) + C", u±(t) = v k {t) = (C, C" G C) 

(cf. [2H]) and that of type B n with 

= U S(*) =^P(*;2wi,2a;2), 

3 

Vk{t) = ^2C j p(t + uj j ;2u; 1 ,2uJ2) - -, (A, C u C £ C) 

introduced by [12] are most fundamental and their integrability and the integrals of higher order 
are established by [23 (23 E2]. Here p(t; 2uj±, 2u>2) is the Weierstrass elliptic function whose 
fundamental periods are 2u>\ and 2lu2 and 

UJq = 0, U>i + W2 + OJ3 = 0. 

Other potential functions are suitable limits of these elliptic potential functions, which is 
shown in [HI [T3"l [33] etc. We will study integrable systems by taking analytic continuations of 
the integrals given in [25j with respect to a suitable parameter, which is done for the invariant 
systems (of type A n -i) by [32] and (of types B n and D n ) by [29] and for the systems of type A n _i 
by [33] • The main purpose of this note is to give the explicit expression of the operators 
P\ , . . . , P n in (jl.2p in this unified way. Namely we construct enough commuting integrals of 
the non- invariant systems from those of the invariant systems given by [241 [29"1 [32] . Such study 
of the systems of types A n _i, B2, B n (n > 3) and D n are explained in Sections El HI El 
respectively. 
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Since the integrals of the system of type A n _\ are much simpler than those of type B n , we 
review the above analytic continuation for the systems of type A n _\ in Section preceding to 
the study for the systems of type B n . There are many series of completely integrable systems 
of type B2, which we review and classify in Section U] with taking account of the above unified 
way. 

We present 8 series of potential functions of type B n in Section [5j There are 3 (elliptic, 
trigonometric or hyperbolic and rational) series of the invariant potentials of type B n whose 
enough integrals are constructed by [25] and [29]. The complete integrability of the remaining 
5 series of the potential functions is shown in Section [5l which is conjectured by [8] (4 series), 
partially proved by [TBI [191 EI] (3 series) and announced by [30] (5 series). The complete 
integrability of two series among them seems to be first established in this note. Note that 
when n > 3, our systems which do not belong to these 8 series of type B n are the Calogero- 
Moser systems with elliptic potentials and the finite Toda lattices of type A^li, whose complete 
integrability is known. 

The main purpose of our previous study in [231 1211 ESI [30] is classification of the completely 
integrable systems associated with classical root systems. In this note we explicitly give integrals 
of all the systems classified in our previous study with reviewing known integrals. 

Since our expression of Pk is natural, we can easily define their classical limits without any 
ambiguity and get completely integrable Hamiltonians of dynamical systems together with their 
enough integrals. This is clarified in Section [71 

In Section[8jwe examine ordinary differential operators which are analogues of the Schrodinger 
operators studied in this note. 



2 Notation and preliminary results 

Let {ei,...,e n } be the natural orthonormal base of the Euclidean space W 1 with the inner 
product 

n 

( x > V) = ^2 x 3 y i for x = (x!,... } x n ), y = (yi, . . . , y n ) G R n . 

3=1 

Here ej = (Sy, . . . , 5 n j) G W 1 with Kronecker's delta 5^. 

Let a G M™\{0}. The reflection w a with respect to a is a linear transformation of W 1 defined 
by w a (x) = x — ^j^x ex for x G W 1 . Furthermore we define a differential operator d a by 



(d a ip)(x) = -rtfix + ta ) 



and then dj = d e . . 



t=o 



The root system S(J3 n ) of type B n is realized in M. n by 

S(^ n _i) + = {ej - &j\ 1 < % < j < n}, 
E(D n )+ = {e i ±e j ; l<i<j<n}, 
£(£„)+ = {e k ; 1 < k < n}, 
E(B n )+ = S( J D n )+ U S(B„)+, 

= {a,-a; a G S(F)+} for F = A n _i, D n or B n . 

The Weyl groups W(B n ) of type B n , W(D n ) of type D n and W{A n _\) of type A n _\ are the 
groups generated by w a for a G T,(B n ), T,(D n ) and S(^4 n _i), respectively. The Weyl group 
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W(A n ^i) is naturally identified with the permutation group & n of the set {1, ... , n} with n 
elements. Let e be the group homomorphism of W{B n ) defined by 



1 if we W(D n ), 
-1 if w £W(B n )\W(D n ). 



e(w) = 

The potential function (|1.3|) is of the form 

ae£(L>„) + aeS(B„)+ 

with functions u a and vp of one variable. For simplicity we will denote 

u a (x) = U- a (x) = u a ((a,x)) for a G 
vp{x) = v-i3(x) = vp({P,x)) for /3 G 

ufj(x) = Ue. ±ej (x), V k (x) = V £k (x). 

Lemma 1. For a bounded open subset U of C, there exists an open neighborhood V of in C 
such that the following statements hold. 

i) The function A sinh" 1 Xz is holomorphically extended to (z,X) G (U \ {0}) x V and the 
function is l/z when X = 0. 

ii) Suppose Re A > 0. Then the functions 

e 2Xt sinh" 2 X(z ± t) and e 4A * (sinh -2 A(z ± t) - cosh -2 A(z ± t)) 

are holomorphically extended to (z,q) G f7 x V with q = e~ 2Xt and the functions are 4e^ 2Xz and 
lQe^ iXz , respectively, when q = 0. 

Proof. The claims are clear from 



A 1 sinhAz = z + ^ (2j + i)! , 

4e" 2A * sinh 2 X(z ±t) = e ±2Xz (l - e~ 2Xt e^ 2Xz ) 2 , 
sinh -2 Xz — cosh~ 2 Xz = 4 sinh -2 2Xz. 



The elliptic functions p and £ of Weierstrass type are defined by 

p(z) = p(z- } a*, a*) = ^ + E ( 71^)2 - ^) . 

C(z) = C(*5 2wi, 2wz) = - + V ( — !— + - + 4), 



where the sum ranges over all non-zero periods 2miu>i + 2m2U)2 (mi, 7712 G Z) of p. The following 
are some elementary properties of these functions (cf. |41j). 

p{z) = p{z + 2uj 1 ) = p(z + 2u2), (2.1) 

C(*) = (2.2) 
(p') 2 = 4p 3 - g 2 p - g 3 = 4(p - ei)(p - e 2 )(p - e 3 ), 

= p(u} u ) for = 1,2,3, W3 = — u>i — and u;q = 0, (2-3) 



Completely Integrable Systems Associated with Classical Root Systems 



5 



1 4 

P( 2 z) = -rS^p{z + u v ) 



u=0 



(12p(z) 2 -g 2 f 
I6p'(z) 2 



(z; 2oj 2 , 2ui) = p(z; 2u>\, 2uj 2 ), 

(ei - e 2 )(ei - e 3 ) 



p(z + wi; 2wi, 2^2) = ei + 



u; 



(z; 2^1,2^2) - ex' 

1.2 



■1A -I 7r, 00) = A 2 sinh" 2 Az + -A 



?(z; 00, 00) = z , 

2w 2 ) = p{z; 2uji,2uj 2 ) + p(z + oji; 2u>i, 2uj 2 ) - e\, 
p(zi) p'(zi) 1 

p{z 2 ) p'{z 2 ) 1=0 if z\ + z 2 + z 3 = 0, 

pfcfc) P'0 3 ) 1 

771 x2 . u _ 2 , ^8nA 2 e" 4nA ^ 
h A smn Az + > „_.,.. cosh 2n\z, 



p(z; 2vi,2u) 2 ) = —— + A 2 sinh" 2 Az + - 
wi , 1 

n 

2/1 00 

C(-i; 2^,2.2) = ^- 



»?1 



n=l 

00 ne~ 4n ^2 



4nAtJ2 I ' 



g = e ™ T = e~ 2Aw2 



n=l 

and A 



7T 



2 A /^To; 1 



Here the sums in (|2.10p converge if 



2Im^> 



UJl \U)\\ 

Let < k < 2m. Then (pTTOj) means 



k 

p z -\ u) 2 ; 2u\,2uj 2 

m 



- * = 4A 2 

UJl 



' 1 00 

--2V- 
12 ^ 1 



nq 



m 

UJl 

2n 



rl k/m p —2\z 
q e 



n n(2-k/m) e 2n\z 



'i_ e -2\ Zq k/ m y ^ 



-,2n 



-,2n 



n=l 



Lemma 2. Let k and m be integers satisfying < k < 2m. Put 

m 



poO; 2wi, 2w 2 ) = p(z; 2wi, 2w 2 ) + 



A 



7T 



2 A /^To; 1 



and i = q 



UJl 

l/m 7riaJ2 /(mui ) 



(2.4) 
(2.5) 

(2.6) 

(2.7) 
(2.8) 

(2.9) 



(2.10) 



Then for any bounded open set U in C x C, there exists a neighborhood V of the origin of C 
such that the following statements hold. 

i) po(z; 2ui, 2uj 2 ) — A 2 sinh" 2 Xz and po(z+ui; 2a; 1, 2^2) + A 2 cosh" 2 Xz are holomorphic func- 
tions of (z, X, q) £ U x V and vanish when q = 0. 

ii) po(z + (k/m)uj 2 ;2u>i,2io 2 ) is holomorphic for (z,X,t) £ U x V and has zeros of order 
min{/c,2m — k} along the hyperplane defined by t = and satisfies 



t po zH a; 2 ; 2a;i, 2a; 2 



in 



k 



AX 2 e~ 2Xz 



(0 < k < m), 



t p [ z-\ uj 2 ;2uji,2u 2 

m 



t=o 



t=o 



8A 2 cosh2Az (k = m), 
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t 



k—2m , 



Z H U>2', 2u>l, 2o>2 

m 



4A 2 e 2Az 



(m < k < 2m). 



t=o 



For our later convenience we list up some limiting formula discussed above. Fix lo\ with 
> and let oj 2 £ M with W2 > 0. Then X = ir/ {2yf— Tuj\) > and 



sinh 2 A(z + wi) = — cosh 2 Xz, cosh 2A(z + u)\) = — cosh 2Az, 



1 



lim A 2 sinh 2 Xz 

A->0 Z z 

lim e 2A| N \ sinh" 2 A(z + iV ) = 4e T2Az , 
lim po(z; 2uj\ , 2^2) = A 2 sinh -2 Xz, 

U)2— >+00 



lim po(^ + wi; 2wi, 2^2) = — A cosh Az, 

1^2^ + 00 



lim e 2rAa;2 po (-2 + rv 2 ; 2wi , 2w 2 ) = 4A 2 e~ 



2Az 



if 



lim e 2AtJ2 p (z + lo 2 ; 2wi, 2w 2 ) = 8A 2 cosh2Az, 
lim e 2 ( 2 - r ) Aw2 p (^ + w 2 ; 2wi,2w2) = 4A 2 e 2Az 

W2— *00 



if 



< r < 1, 



1 < r < 2. 



3 Type A n _! (n > 3) 

The completely integrable Schrodinger operator of type A n _\ is of the form 



I d 2 

p = — 2 Y.q^+ E 

j=l 3 \<i<j<n 



Xi Xj j 



Denoting 



u ei-ej{ x ) — U e .- ei (x) — U-ixi Xj), 



we put 

p k = E 



E U w{e 1 - 



2 u u\(k-2u)\(n-k)\ 

E E U *i«2 ' ' ' "ta^-ita^^+i " ' ^ 
0<i/<fc/2 



1 k 



(2.11) 
(2.12) 
(2.13) 
(2.14) 
(2.15) 
(2.16) 
(2.17) 
(2.18) 



(3.1) 



according to the integrals given in \25\ I32|. We will examine the functions u-(t) which satisfy 
[Pi,Pj] =0 for 1 < i < j < n. (3.2) 

Here we note that 

P = P 2 -\pI Pi = di + --- + <9 n , 
P2 = E did i + E u ij( x i- x i)> 

l<i<j<n 1 <i<j? <n 

n 

P 3 = didjd k + J2 E u ij( x i- x j) d k- 



l<i<j<k<n 



k=l l<*<J<n 
i^k, j^k 
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Pa= Y d i d Ade+ Y S u iA d e+ Y ( u ij u ke + u ik u je + u ie u jk) 

l<i<j<k<£<n l<k<£<n l<i<j<n l<i<j<k<£<n 

= Y d i d A & £ + Y U ij d k & i + Y U ij U kV 

p 5 = Y didjdkdtd m + Y u ijdkdid m + Y u ij u k£ d m, 

p 6 = Y didjd k d £ d m d u + Y u~-d k d(d m d v + Y u ij u ki^m^ + Y u ij u k£ u mw 



Since W{A n _\) is naturally isomorphic to the permutation group & n of the set {1, . . . , n}, we 
will identify them. In |25l I32j . the integrable potentials of type A n ^i which are invariant under 
the action of S n are determined and moreover (|3.2p with (|3.ip is proved. They are 

u ei - ej (x) = u((ei - ej,x)) (3.3) 

with an even function u and 

(Ellip-A„_i) Elliptic potential of type A n _i: 

u(t) = Cp (t; 2wi,2a;2), 

RE(A n -r,xx, . . . ,x n ;C,2ui,2u)2) = C Y Po(°°i ~ ^jS 2 ^i, ^2), 

l<i<j<n 

(Trig-A„_i) Trigonometric potential of type A n _i: 

u(t) = C sinh -2 Xt, i?y(A n _i; x±, . . . , x n ; C, A) = C Y2 sinh -2 A(a 

l<i<j<n 



(Rat-A^-i) Rational potential of type A. 



n-l- 



C Q 
u(t) = ~2, R R (A n -x;xi,...,x n ;C) = Y 



l<i<j<n v 4 J ' 

We review how the integrability of (Ellip-^n-i) implies the integrability of other systems. 
Since it follows from (|2.14p that 

lim R E ( A n _i;x; ^,2^1, 2^) = R T (A n -x; x; C, A), 



UJ2— >00 



c 

u(t) = lim — ^po{t; 2lu±, 2ui2) = C sinh -2 Xt, 

the integrability (|3.2p for (Trig-^4 n _i) follows from that for (Ellip-^4 n _i) by the analytic conti- 
nuation of u ei _ ej (x) and P k with respect to q (cf. (|2.10p . (|3.ip . f|3.3|) and Lemma[2]i)). 
The integrability for (Rat-^4 n _i) is similarly follows in view of Lemma [1] with 



lim R T (A n _i;x; X 2 C, X) = R R (A n ^i; x; C), 
A — >0 

u(t) = lim A 2 Csinh" 2 Ai. 



This argument using the analytic continuation for the proof of (|3.2h is given in 

As is shown |33j . there are two other integrable potentials of type A n _i to which this argument 
can be applied. 
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(Toda-^ 1 2 1 ) Toda potential of type A^^. 



n-l 

.(1) . 



Rl{A^ ] _i, z; C, A) = Ce x{x *- x > +l) + Ce x ^~ Xl \ 



i=l 



(Toda-74„_i) Toda potential of type A n . 



R L {A n _ x - x; C,X) = J2 Ce x ^~ x ^. 



i=l 

The integrability (|3.2p for these potentials similarly follows in view of Lemmas HJ [2] and 

/ . 2uJ 2 2kuJ 2 e (4/n)Aa. 2 

hm R E Ai-i; x x , . . . , x k , ...,x n - 2lo 2 ; — rrn — C, 2lo±, 2uj 2 

uj 2 ^oo \ n n AA Z 

= R L (A^_ l ;x;C,-2X), 

3 (4/n)Aoj 2 



u et -e (x)= lim — — - — Cpolxi- xj -\ ;2ui,2u 2 

{ Ce- 2X{ - Xi ~ Xi +^ if 1< j = % + 1 < n, 

Ce -2\(x n - Xl ) if i = l and ■ = ^ 

if 1 < % < j < n and j — i ^ 1, n — 1 

and 

lim i?r A n _i; x± — N, . . . ,x n — nN; — - — C, A = i? L (^n-i! C, -2A), 
./V— >oo y 4 y 

e 2AJV 

u 6i _ e (a;) = lim — — C sinh -2 A(xi - xj + (j - i)N) 

r c . e -2A(x i -x i+1 ) j£ i< j = j + 1 < n> 

1 if 1 < i < j < n and j ^ i + 

respectively, if Re A > 0. The restriction Re A > is removed also by the analytic continuation. 

Thus the following theorem is obtained by the analytic continuation of the integrals (|3.1|) of 
(Ellip-^4 n _i) whose commutativity ()3.2[) is assured by [52"] . 

Theorem 1 [101 llll 128] I32 |. [33 ] 134] . etc.). TTie Schrodinger operators with the potential 

functions (Ellip-A„„i), (Trig-j4 n _i), (RaWl n _i), (Toda-^ 1 2 1 ) and (Toda-^4 n _i) are completely 
integrable and their integrals are given by (|3.ip with u ei - ej (x) in the above. 

Remark 1. i) There are quite many papers studying these Schrodinger operators of type A n _\. 
The proof of this theorem using analytic continuation is explained in [55] . 

ii) The complete integrability (|3.2p for (Ellip-A n „ 1 ) is first established by [32] Theorem 5.2], 
whose proof is as follows. The equations = [P2)-Pfc] = for k = 1, ...,n are easily 
obtained by direct calculations with the formula (|2.9p (cf. [321 Lemma 5.1]). Then the relation 
[P 2 , [Pi,Pj]] = and periodicity and symmetry of Pk imply [Pi,Pj] = (cf. [52 Lemma 3.5]). 
Note that the proof of the integrability given in |28[ § 5 Proposition 2 and Appendix E] is not 
correct as is clarified in [32] Remark 3.7] (cf. [33] § 4.2]). 

Note that the complete integrability for (Trig-j4 n _i) is shown in 

iii) If Re A > 0, we also have 

lim R L (A® l]Xk + kN;Ce- 2XN -,-2\) = R L {K-i\ »5 C, -2A). 
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iv) Note that 

R L (A n ^; x k + N k ; C, A) = ^ C e x ^~ N ^ ■ e x ^~ x ^. 



n-l 



i=l 

Hence e xlyXl ~ X2 ^ — e x ^ X2 ~ X3 ^ gives the potential function of a completely integrable system of 
type A n _i with n = 3 but the potential function 

lim \-i( e Hxi-*2) _ e K*2-x a )\ _ 2x2 + X3 

A->0 v 7 

does not give such a system because it does not satisfy (|9.3p in Remark [T31 

Considering the limit of the parameters of the integrable potential function, we should take 

care of the limit of integrals. 

v) Let P n (t) denote the differential operator P n in (|3.ip defined by replacing u~- by viZ = u ij+t 

with a constant i £ C. Then 

p «(*)= 2 wr p "- 2fctfc with Po = 1 ' 

0<fc<n/2 

[J o n(«),i'n(<)]=0 for s,t€C. (3.4) 
In fact, the term li^itjj^ ■ ■ ■ u 2j-i 2j®'2j+id2j+2 • • • d n -2k appears only in the coefficient of t k in 



the right hand side of (|3.4p and it is contained in the terms 

U7 i ■■■U~ i U7n ■ ■ ■ Uoa 1 o^do-i+i ■ ■ ■ <9n-2fc 

«n-2fc + l«n-2fc + 2 «n-l«n 12 2j-l,2j 47 T 1 n ZK 

of P n (t), where the number of the possibilities of these Ui n _ 2k+1 i n _ 2k+2 ■ ■ -u7 li is (2ft)!/(2 fc ft!) 
because {i n -2k+i, «n-2fc+2, • ■ • , vl = {n - 2ft + 1, . . . , ra}. 
vi) Since 

= (n - ft + H \- x n ] for ft = 2, n, 

[Pfc, P2] = implies [Pk-i, P2] = by the Jacobi identity. Here we note that 

[ U 12 U M • • • U 2j-l,2j 9 2j+l • • • d k ~ld„, X V ] = Ui 2 U3 4 • • • U^-l^j^Zj+l ' ' ' ^fc-l 
for ^ = ft, ft + 1, . . . , n. 

vi) The potential functions of (Trig-^4 n _i), (Rat-^4 n _i) and (Toda-^4 n _i) are specializations 
of more general integrable potential functions of type B n (cf. Definition [5]). 

In the following diagram we show the relations among integrable potentials of type A n ^i by 
taking limits. 

Hierarchy of Integrable Potentials of Type A n -i (n > 3) 

Toda-A^ 1 ^ — ► Toda-A n _i 

/ / 
Ellip-A n _i Trig-A n „! Rat-A n _! 
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4 Type B 2 

In this section we study the following commuting differential operators P and P 2 - 
1 {_&_ 3_ 2 



+ S with ord S < 4, 



dx 2 dy 2 

[P,P 2 }=0. (4.1) 

The Schrodinger operators P of type P>2 in this section are known to be completely integrable. 
They are listed in [23| 12^1 130] . We review them and give the explicit expression of P2 ■ 

First we review the arguments given in [231 124] . Since P is self-adjoint, we may assume P2 is 
also self-adjoint by replacing P2 by its self-adjoint part if necessary. Here for 

A = ^a lJ (x,y)Q^J 
we define 

*A = Yi-iy^^^aaix^) 
^ ' dx l dyi n ,yj 

and A is called self-adjoint if A = A. 

Lemma 3 ([23]). Suppose P and P2 are self-adjoint operators satisfying ([4. II) . Then 
R(x, y) = u + (x + y) + u~(x - y) + v(x) + w(y), 

(q2 \ ^ q2 q2 

Q^+ u ~( x -y)- u+ ( x +y)j ~ 2w{y)-^-2v{x) — +Av{x)w{y)+T{x,y) (4.2) 

and the function T(x, y) satisfies 

— ! — = 2(u + {x + y)-u~(x- y)) — h Aw{y) — {u + (x + y)-u~(x- y)), 

ox oy ay 

q = 2(u + (x + y)-vT{x- v))-£r + 4t;(x) — (n+(x + y)-u~(x- y)). (4.3) 

Conversely, if a function T(x,y) satisfies (|4.3p for suitable functions (t) , v(t) andw(t), then 
(|4.1|) is valid for R(x,y) and P2 defined by (|4.2|) . 

Remark 2. i) If w(y) = 0, then T(x,y) does not depend on x. 

ii) The self-adjointness of P2 and the vanishing of the third order term of [P, P2] imply that 
P2 should be of the form (|4.2p with a suitable function T(x,y). Then the vanishing of the first 
order term implies (|4.3|) . The last claim in Lemma [3] is obtained by direct calculation. 

Since T(x, y) satisfying (|4.3j) is determined by (u~ ,u + ;v, w) up to the difference of constants, 
we will write T{u~ , u + ; v , w) for the corresponding T{x,y) which is an element of the space of 
meromorphic functions of (x,y) modulo constant functions and define Q(u~ ,u + ;v,w) by 

T(u~ ,u + ;v, w) = 2(u~(x — y) + u + {x + y))(v(x) + w(y)) — AQ(u~ ,u + ;v, w). (4-4) 

Note that T(u~ ,u + ;v,w) and Q(u~ , u + ; v , w) are defined only if the function T(x,y) satis- 
fying (|4.3|) exists. The following lemma is a direct consequence of (|4.3|) and this definition 
of Q. 
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Lemma 4 ([23]). Suppose T(u ,u + ;v,w) and T{u i ,uf]Vj,Wj) are defined. Then for any C, 
Ci, Cj G C the left hand sides of the following identities are also defined and 

T{u~{t) + C, u + (t) + C; v(t), w(t))) = T(u'(t),u + (t); v(t),w(t)), 
Q(u~ {t),u + (t);v(t) + C,w(t) + C)) = Q(u-(t),u + {t);v{t),w{t)), 
Q(u'(Ct),u + (Ct); v(Ct),w(Ct)) = Q(u~ (t),u + (t); v(t),w(t))\ x ^ C x, y^Cy, 

(2 2 2 2 \22 

Y AiU i> Y AiU t' Y C i v ii Y C i w 3 )=YY A i°jQ( u i ' u ti v i' w 3)- 
i=l i=l j=l j=l ) i=l j=l 

Hence the left hand sides give pairs P and P2 with [P, P2] = 0. 
For simplicity we will use the notation 

Q(u~ , u + ; v) = Q{u~ , u + ; v , v ), Q(u;v,w)=Q(u,u;v,w), Q(u;v) = Q(u,u;v,v). (4.5) 

The same convention will be also used for T{u~ ,u + ;v,w). The integrable potentials of 
type P>2 in this note are classified into three kinds. The potentials of the first kind are the 
unified integrable potentials which are in the same form as those of type B n with n > 3, which 
we call normal integrable potentials of type B<i- 

The integrable potentials of type B2 admit a special transformation called dual which does 
not exist in B n with n > 3. Hence there are normal potentials and their dual in the invariant 
integrable potentials of type B2. Because of this duality, there exist another kind of invariant 
integrable potentials of type B2, which we call special integrable potentials of type B2. 

In this section we present (R(x,y),T(x,y)^ as suitable limits of elliptic functions as in the 
previous section since it helps to study the potentials of type B n in Section [5j We reduce the 
complete integr ability of the limits to that of a systems with elliptic potentials. But we can also 
check (|4.3p by direct calculations (cf. Remark [5]). 

4.1 Normal case 

In this subsection we study the integrable systems (]4.ip with (|4.2p which have natural extension 
to type B n for n > 3 and have the form 

3 3 
u~(t) = Auo(t), u+{t) = Au+{t), v(t) = C jVj (t), w(t) = C jWj {t) 

i=o i=o 

with any A, Co, C±, C2, C3 G C. These systems are expressed by the symbol 
(( n o )' ( n o~); ( v o,v 1 ,v 2 ,v 3 ), (w ,wi,w 2 ,w 3 )). 

The most general system is the following (Ellip-i?2) defined by elliptic functions, which is called 
Inozemtsev model [12] , 

Theorem 2 (-B2: Normal case, \12 \ \24 \ [23 | I30| etc.). The operators P and P2 defined by the 
following pairs of R(x,y) and T(x,y) satisfy (|4.1|) and (|4.2[) . 

(Ellip-£ 2 ): ((p{t)}; (p(t), p{t + p(t + u 2 ), pit + w 3 )» 
3 3 

j=0 j=0 
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u (x-y) = Ap(x-y), u + (x + y) = Ap(x + y), 

3 

R(x,y) = A(p(x -y) + p(x + y)) + s ^C j {p(x + uj) + p(y + Uj)), 



T(x, y) = 2A(p(x -y) + p(x + y)) ^ Cj (p(x + Uj) + p{y + Uj)j 

3 

- 4A^2Cjp(x + ojj)p(y + cjj). 

j=0 

(Trig-5 2 ): ((sinh -2 At); (sinh -2 At, sinh -2 2At, sinh 2 At, sinh 2 2At)) 

v(x) = Cq sinh" 2 Xx + C\ cosh" 2 Ax + C 2 sinh 2 Ax + -C3 sinh 2 2Ax, 

w(y) = Co sinh" 2 Ay + C\ cosh" 2 Xy + C2 sinh 2 Xy + -C3 sinh 2 2Ay, 

u"(x — y) = A sinh" 2 A(x — y), u + (x + y) = A sinh" 2 A(x + y), 
R(x, y) = ^(sinh -2 A(x — y) + sinh" 2 A(x + y)) + Co (sinh -2 Ax + sinh" 2 Ay) 
+ C\ (cosh" 2 Ax + cosh -2 Ay) + C2 (sinh 2 Ax + sinh 2 Ay) 



^3(S 

-2 \/„, „.\ , „:„u-2 



+ ^C 3 (sinh 2 2Ax + sinh 2 2Ay) . 



T(x, y) = 2^4(sinh A(x — y) + sinh A(x + y)) 

x ^Co(sinh" 2 Ax + sinh" 2 Ay) + Ci(cosh" 2 Ax + cosh" 2 Ay) 

+ C 2 (sinh 2 Ax + sinh 2 Ay) + ^C 3 (sinh 2 2 Ax + sinh 2 2 Ay) 

— 4A^Co sinh -2 Ax • sinh" 2 Ay — C\ cosh" 2 Ax • cosh" 2 Ay 
+ C3(sinh 2 Ax + sinh 2 Ay + 2 sinh 2 Ax • sinh 2 Ay) 

(Rat-£ 2 ): ((t -2 ); (t -2 , t 2 , t 4 , t 6 )) 

v{x) = Cox- 2 + C x x 2 + C 2 x 4 + C3X 6 , w(y) = C y- 2 + C lV 2 + C 2 y 4 + C 3 y 6 , 
iT(x-y) = - — ^—2, u + (x + y) 



(x-y) 2 ' (x + y) 2 ' 

= J^Zyf + J^yJ +Co(x~ 2 + y- 2 )+C 1 {x 2 +y 2 ) + C 2 {x A +y 4 ) + C 3 (x 6 +y 6 ) 

T(x, y) = 2(u - (x - y) + u + (x + y))(v(x) + w(y)) - 4A( + C 2 (x 2 + y 2 ) 

2„,2 1 r< J2„,2(J2 1 „,2\ 1 op „,4„,4 

+ C 3 (x 4 + y 4 + 3x 2 y 2 ) ) =8A-' 



2 2^ \ _ „ 2Co + 2CixV + C 2 x 2 y 2 (x 2 + y 2 ) + 2C 3 x 4 y 4 



(x 2 — y 2 ) 2 

(Toda-D^-bry): ((cosh 2At); (sinh -2 At, sinh -2 2At), (sinh -2 At, sinh -2 2At)) 

v(x) = Co sinh -2 Ax + C\ sinh -2 2Ax, w(y) = C 2 sinh -2 Ay + C3 sinh -2 2Ay 
u~(x — y) = ^4cosh2A(x — y), u + (x + y) = Acosh2A(x + y), 
R(x, y) = A cosh 2A(x — y) + A cosh 2A(x + y) 

+ Co sinh -2 Ax + C\ sinh -2 2Ax + C 2 sinh -2 Ay + C 3 sinh -2 2Ay, 
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T(x,y) = 8A(C cosh 2Xy + C 2 cosh 2\x). 

(Toda-S^-bry): ((e~ 2At ); (e 2A * , e 4A <), (sinh~ 2 At, sinh" 2 2At)) 

= C e 2A:c + Cie 4Ax , w(y) = C 2 sinh" 2 Ay + C 3 sinh~ 2 2Ay, 
u - (x-y) = Ae~ 2X{ - x - y) , u + {x + y)= Ae~ 2X(x+y) , 

R{x, y) = Ae- 2X{x - y) + Ae~ 2X{x+y) + C e 2Xx + de iXx + C 2 sinh~ 2 Ay + C 3 sinh" 2 2Ay, 
T(x, y) = AA(C cosh 2Ay + 2C 2 e~ 2Xx ) . 

(THg-4 1} -bry): ((sinh^ 2 At), 0; (e~ 2A * , e" 4At , e 2A4 , e 4A *» 

v{x) = C e~ 2Xx + de- iXx + C 2 e 2Xx + C 3 e 4Xx , 
to(y) = C e~ 2Xy + C ie - AX y + C 2 e 2Ay + C 3 e AXy , 
u~(x — y) = ^4sinh~ 2 X(x — y), u + (x + y) = 0, 

= A sinh" 2 \(x - y) + Co{e~ 2Xx + e~ 2Xy ) + d(e- 4Xx + e~ iXy ) 

+ C 2 (e 2Xx + e 2Xy ) + C 3 (e AXx + e AXy ), 
T(x, y) = 2A sinh" 2 X(x - y)(C (e- 2Xx + e~ 2Xy ) 

+ 2C 1 e~ 2X( - x+y *> + C 2 (e 2Xx + e 2Xy ) + 2C 3 e 2X{x+y) ). 

(Toda-Cf ] ) : ((e~ 2A * ) , 0; (e 2A * , e 4At ) , (e~ 2A * , e~ 4Ai ) ) 

v ( x ) = C e 2Xx + Cf Xx , w(y) = C 2 e~ 2Xy + C 3 e~ 4Xy , 
u - (x-y) = Ae~ 2X{ - x - y) , u + (x + y) = 0, 
R(x, y) = Ae~ 2X{x - y) + C e 2Xx + Cf Xx + C 2 e~ 2Xy + C 3 e~ iXy , 
T(x,y) = 2A(C e 2Xy + C 2 e~ 2Xx ). 

(Rat-^i-bry): ((i~ 2 >, 0; (t, t 2 , t 3 , t A )) 

v{x) = C x + C lX 2 + C 2 x 3 + C 3 x\ w(y) = C y + C\y 2 + C 2 y 3 + C 3 y\ 

u~(x-y) = - f — ^-Tn, u + (x + y) = 0, 
(x - yY 



R(x, y) = - +C (x + y)+ d(x 2 + y 2 ) + C 2 (x 3 + y 3 ) + C 3 (x 4 + y^ 
(x - y) z 

2A 



T(x,y) = -^—(C (x + y)+C 1 (x 2 + y 2 )+C 2 xy(x + y) + 2C 3 x 2 y 2 ). 
[x - yY 

Remark 3. For example, ((t _2 ),0; (t, t 2 , t 3 , t 4 )) in the above (Rat-^-bry) means 

u~(t) = At" 2 , u + (t) = 0, v(t) = w{t) = C t + Cit 2 + C 2 t 3 + C 3 t 4 

with a convention similar to that in (|4,5D . 

We will review the proof of the above theorem after certain remarks. 

Remark 4. All the invariant integrable potentials of type B 2 together with P 2 are determined 
by [24^125] . They are classified into three cases. In the normal case they are (Ellip-^), (Trig- 
and (Rat-£?2) which have the following unified expression of the invariant potentials given by 
Lemma 7.3], where the periods 2oj\ and 2u;2 may be infinite (cf. (|2.6p and (|2.7p ) 

R( X , y) = A P(I -y) + A P ( x + y) + + C ^ ± ™*? ± CK*) ± C, 
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+ 



C A p(y) 4 + C 3 p(y) 3 + C 2 p(y) 2 + C lP (y) + C 



P'{yf 



c 3 



T(x,y) = 4A(p(x) - p(y))- 2 [C 4 p(xVp(yy + -fp{xYp{y) 

c c c 

+ -Yp{x)p{y? + C 2 p{x)p{y) + —p{x) + —p{y) + C 

This is the original form we found in the classification of the invariant integrable systems of 
type B n (cf. [25]). Later we knew Inozemtzev model and in fact, when the periods are finite, 
(I2.3h and (|2.5p show that the above potential function corresponds to (Ellip-i^)- 

y 2 



When uji = ll>2 = oo, p(t) = t 2 and (p(x) — p{y)) 2 = x A y 4 (x 2 
R(x,y) = - — + 



2 and 



A A I 

7 ttt + t ttt + t(C 4 2;- 2 + C 3 + C 2 x 2 + Cix 4 + C x 6 ) 

[x — y) z [x + y) z 4 



C 4 y" 2 + C 3 + C 2 y 2 + C lV A + C y 6 ), 

2 „.2\-2/r,n , n (Ji , „.2\ , o/-v „,2„,2 , „,2„.2/ 2 , „,2\ , o/^ „,4„.4n 



T(x,y) = 2,4(x^ - 2/T (2^4 + C 3 (x^ + y z ) + 2C 2 xV + CixV(^ + y z ) + 20,* V)- 
We review these invariant cases discussed in |24[ [32] . Owing to the identity 



u 


+ )v' + iv((u~ 


r 


-(u 


^)')+d y {2{u- 




v v' 


i 




v —v' 1 


2 


w —w' 


i 


+ 2 


10 —w' 1 




vT -{u-y 


i 




u+ (u+y i 



and (|2.9p . the right hand side of the above is zero and we have (|4.3j) when 

u - = Cp(x-y) + C, u + = Cp(x+y) + C , v = Cp{x)+C and w = Cp(y) + C 

with T(x,y) = 2(n™ + u Jr )(v + w) — Avw. Hence with Q(p(t); p(t)) = p(x)p(y), the function 
T(x,y) given by (j4.4H satisfies (14. 3p with the above u^, v and Using the transformations 
(x, y) (x + ujj,y + ujj), we have 



(x + Uj)p(y + Wj) = Q(p(t + uj); p(t + Wj)) = Q(p(t); p(t + 



(4.6) 



for j = 0,1,2,3 because the function p(x ± y) does not change under these transformations 
(cf. (|2.ip ). Thus we have Theorem [2] for (Ellip-i? 2 ) in virtue of Lemma [3] and Lemma 01 
Here we note that p may be replaced by po. 

By the limit under uj 2 — > oo, we have the following (Trig-i? 2 ) from (Ellip-i? 2 ). See the proof 
of [291 Proposition 6.1] for the precise argument. 



(Trig-5 2 ): 

C;(sinh -2 At; sinh -2 At) = sinh -2 Xx ■ sinh -2 Ay, 
<5(sinh -2 Xt; cosh -2 Xt) = — cosh -2 Ax • cosh -2 Ay, 
Q(sinh -2 Xt; sinh 2 At) = 0, 

sinh - At; — sinh 2At I = sinh Ax + sinh Xy + 2 sinh Ax • sinh Xy. 



(4.7) 
(4.8) 
(4.9) 

(4.10) 



The equations (|4T7| . (gSD and ([43]) correspond to (gTll}, P7T3]) and §1~TI\ , respectively. More- 
over (|2.8p should be noted and (|4.10p corresponds to (|2.17p with replacing (u>i, X) by (cji/2, 2A). 
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By the limit under A — > 0, we have the following (Rat-i?2) from (Trig-i^) as was shown in 
the proof of |29t Proposition 6.3]. Here we note (|2.13|) and 

cosh -2 At • sinh 2 At = 1 — cosh" 2 At, 

cosh -2 At • sinh 4 At = — 1 + cosh -2 At + sinh 2 At, 

cosh -2 At • sinh 6 At = 1 - cosh -2 At - 2 sinh 2 At + - sinh 2 2At, 

4 

lim X~ 2j cosh -2 At • sinh 2j ' At = t 2j for j = 1, 2, 3, 

A^O 

1 1 _2(x 2 + y 2 ) 



(x — y) 2 ' (x + y) 2 (x 2 — y 2 ) 2 
The result is as follows. 

(Rat-£ 2 ): 

Q(t~ 2 ;t~ 2 ) = x -2 y -2 , 

T(t -2 ; t -2 ) = 4(x -2 + y -2 )((x - y) -2 + (x + y) -2 ) - 4x -2 y -2 
4(x 2 + y 2 ) 2 - 4(x 2 - y 2 ) 2 _ 16 

x 2 y 2 {x 2 — y 2 ) 2 (x 2 — y 2 ) 2 ' 

0, 

4(x 2 + y 2 ) 2 16x 2 y 2 



Q(t~ 2 : 


,t 2 ) 


T(t -2 : 


t 2 ) 


Q(t~ 2 : 


A 


T(t -2 , 




Q(t~ 2 : 




T(t- 2 , 


t 6 ) 



(x 2 — y 2 ) 2 (x 2 — y 2 ) 2 
Ji , „.2 



+ 4, 



4(x 2 +y 2 )(x 4 + y 4 ) 2 2 , _8xV(x 2 +y 2 ) 

<i[x + y 



(x 2 — y 2 ) 2 (x 2 — y 2 ) 2 

4 i 4 i o 2 2 

x + y + 3x y , 

4^ + ^ + ^ 4(i1 + 3iV+!)1) _ 16xV 



(x 2 — y 2 ) 2 (x 2 — y 2 ) 2 

This expression of T(x,y) for (Rat-.E>2) is also given in Remark |H Note that we ignore the 
difference of constants for Q and T. 

Proof of Theorem [2], The three cases (Ellip-i?2), (Trig-i?2) and (Rat-i?2) have been explai- 
ned. Note that if u , v, w and T(u~ ,u + ;v,w) (or Q(u~, u + ; v, w)) are defined and they have 
an analytic parameter, Lemma [3] assures that their analytic continuations also define P and P2 
satisfying [P, P 2 ] = 0. 

(Toda-D^'bry) <— (Ellip-£>2): Replacing (x,y) by (x + a;2,y), we have 

( e 2Aw 2 1 1 

— r p {t + uj 2 ); t^Po(* + w 2 ), -r^poW 

lim ""ST^PoO + oj 2 )po(y) = cosh2Ax • sinh -2 Ay, 



u>2— *oo 



8A 4 



/ g2Aw 2 I I 

Q(cosh2At;0,cosh -2 At) = lim Q — -^-poO+u^); --To Po(*+^i + uj 2 ); --— po(t+wi 

o; 2 -»oo y 8A Z A z A2 

e 2Aw2 

= ~ lim -sri-Po^ + w i + ^2) • po(y + ^l) 
oj 2 ->oo 8A 4 

= — cosh 2 Ax • cosh -2 Ay, 
Q(cosh 2At; sinh -2 At, 0) = sinh -2 Ax • cosh 2 Ay, 
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Q(cosh 2At; cosh 2 At, 0) = — cosh 2 Ax • sinh 2 Ay. 
Hence 

T(cosh 2At; 0, sinh" 2 At) = 2(cosh 2A(x + y) + cosh 2A(x - y)) ■ sinh" 2 Xy 

— 4 cosh 2Xx ■ sinh" 2 Xy = 8 cosh 2Xx, 

T(cosh 2Xt; 0, cosh" 2 At) = 2 (cosh 2A(x + y) + cosh 2A(x — y)) • cosh -2 Ay 

+ 4 cosh 2Ax • cosh -2 Ay = 8 cosh 2Ax, 

T(cosh 2At; 0, sinh -2 2Xt) = T(cosh 2At; sinh -2 2At, 0) = 0, 

T(cosh 2At; sinh -2 At, 0) = 8 cosh 2Ay. 

(Toda-S^-bry) <— (Toda-D^-bry): Replacing (x,y) by (x — N,y), we have 



T(e -2A * ; 0, sinh -2 At) = lim T{2e~ 2XN cosh 2A(t - N)\ 0, sinh -2 At) 

TV^oo 

-2ATV „„„, o\/„ Art o„— 2Aa 



= lim 16e"^ v cosh2A(x - TV) 

TV^oo 

T(e -2A *;0,sinh -2 2At) = lim T(2e -2AAr cosh 2A(t - iV); 0, sinh -2 2At) = 0, 

TV— >oo 

T(e -2A *;e 2A ',0) = lim T[2e~ 2XN cosh 2X(t-N); -e 2XN sinh -2 A(t-iV), 0) = 4cosh 2Ay, 

TV— >oo y 4 y 

T(e -2A *;e 4A *,0) = Jim T^2e -2AAf cosh 2A(t - iV); i e 4A7V sinh -2 2A(t - N), =0. 

(Trig-C^) <- (Trig-S^-bry): Replacing (x,y) by (x + N,y + N), 
T(e- 2Xt ,0;e 2Xt ,0)= lim T(e~ 2Xt , e - 2 M*+2/v). e -2ATV e 2A(t+TV) ? Q) 

TV— >oo 

= lim e -2AJV cosh 2A(y + tV) = 2e 2Ay , 

TV— >oo 

T(e -2A ',0;e 4At ,0) = lim T(e -2A ', e -2A ( t+27V ); e -4A7V e 4A (' +JV ), 0) = 0. 

TV— >oo 

By the transformation (x, y, A) i— > (y, x, —A), we have 

T(e -2At , 0; 0, e -2At ) = 2e~ 2Xx , T{e~ 2Xt , 0; 0, e -4A *) = 0. 

(Trig-.A^-bry) <— (Trig-i?2): Replacing (x,y) by (x + tV, y + iV), 

Q(sinh -2 At,0;e 2A ') = lim Q (sinh -2 At, sinh -2 A(t + 2iV); -e 2XN sinh -2 A(t + NU 
= lim \e 2XN sinh -2 A(x + iV) sinh -2 A(y + N) = 0, 

TV— >oo 4 



T(sinh -2 At, 0; e 2A ') = 2(e 2Al + e 2Xy ) sinh -2 A(x - y), 



Q(sinh -2 At, 0; e 4A *) lim Q(sinh -2 At, sinh -2 A(t + 2tV); 4e -4A7V sinh 2 2A(t + N)) 

TV— >oo 

-4ATV/ ■ i 2 w , AT \ , „;„u2 



16 lim e ~ 4A7V (sinh A(x + N) + sinh A(y + iV) 

TV^oo 



+ 2 sinh 2 A(x + N) sinh 2 A(y + N)) = 2e 2X{x+y) , 
T(sinh -2 At, 0; e 4A *) = 2 sinh -2 A(x - y)(e 4Xx + e 4Xy ) - 8e 2X( - x+y ^ 

= 2sinh -2 A(x - y)(e 4Ax + e 4Xy - e ^+y) ( e K*-y) _ e -A(*-w))2) 
= 4e 2A ( x+ ^sinh -2 A(x-y). 
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(Rat-Ai-bry) <— (Trig-^4i-bry): Taking the limit A — ► 0, 

Q(t~ 2 ,0;t) = lim Q ( X 2 sinh" 2 At, 0; ^r(e 2Xt - 1) ) = 0, 

T(t- 2 ,0;t) = ^±|, 
[x - yy 

Q(t~ 2 ,0;t 2 ) = limqf A 2 sinh- 2 At,0;4^(e 2A ' + e~ 2Xt - 2) ) = 0, 
a^o \ AX 2 

T(t- 2 ,0;t 2 ) = 2- x2 + y2 



(x - y) 2 ' 

Q(r 2 , 0; t 3 ) = lim Q ( A 2 sinh" 2 At, 0; -^(e 4Af " 3e 2A * - e~ 2A * + 3) 
a->o \ SA^ 

= }2iJ^ 2i(x+ " , - 1 > = ^ + A 

T(^0; ( ») = 2fi±4-2(, + !/ )=2^±|, 
(x - y) 2 (x - y) 2 

Q(t~ 2 , 0; t 4 ) = lim Q ( A 2 sinh" 2 At, 0; -Ax(e 4At + e" 4A * - 4e 2Xt - 4e~ 2Xt + 6) 
A-^o \ 16A 4 

= ft T^ {e2X[x+v) + e " 2A(x+y) - 2) = ^ + y)2 ' 

r(t- 2 ,0;t 4 ) = 2^±^-2(x + y) 2 - 4 ^ V 



(x-y) 2 {x-y) 2 ' 
Thus we have completed the proof of Theorem [2] by using Lemma [3] and Lemma 2J ■ 

Remark 5. Theorem [2] can be checked by direct calculations. For example, Remark [2] and the 
equations 

2(ee~ 2X( - X+v) - e -2A(s-i/)w e 2Axy + , ee -2\(x+y) _ e -2X(x-y)\ 

dx 

= 4X(ee- 2Xy - e 2Xy ) - 8X{ee~ 2Xy - e 2Xy ) = ^-(2(ee- 2Xy + e 2Xy )), 

dy 

2 , £e -2X( X+ y) _ e -2X(x-y)^ e AXxy + 4g 4Ax 9 ,-2X(x+y) _ e ~2X(x~ y)) = q 

ax 

with e = 1 give T(e~ 2Xt ; e 2Xt , 0) and T(e~ 2Xt ; e 4Xt , 0) for (Trig-5 2 -bry). Moreover the functions 
T(e~ 2A ',0;e 2A ',0) and T(e~ 2Xt , 0; e 4A * , 0) for (Toda-C^) also follow from these equations with 
e = 0. 

4.2 Special case 

In this subsection we study the integrable systems (|4.ip with (|4.2p which are of the form 
i?(x, y) = u~(x — y) + u + (x + y) + t>(x) + u>(y), 



7/ 



i i 

"(*) = E <>^)' = E 

i=o j=o 

with Ao, Ai, Co, C\ G C. The most general system (Ellip-^-S) in the following theorem is 
presented by [23] as an elliptic generalization of (Trig-I?2-S) found by [32] . 
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Theorem 3 (B 2 ■ Special case, [23J [Ml ED] etc.). The operators P and P 2 defined by the 
following pairs of R(x,y) andT(x,y) satisfy (|4.1|) and (|4.2|) , 

(Ellip-5 2 -S): ((p(t',2ui,2u 2 ), p(t;ui,2u 2 )); (p{t;u\, 2u 2 ), p(t; u±, uj 2 ))) 

v(x) = Cop{x;u>i,2u2) + C\p{x;u\, u 2 ), w(y) = C p(y; 2u 2 ) + Cip(y; ujx, uj 2 ), 
u~(x -y) = A p(x - y; 2ut,2u 2 ) + A 1 p(x - y;ui,2u 2 ), 
u + (x + y) = A Q p{x + y;2ui,2u 2 ) + A\p{x + y; wi, 2u> 2 ), 
R(x,y) = A p(x - y; 2ut,2u; 2 ) + A p(x + y; 2ui, 2u 2 ) 
+ A\p{x - y; ui,2u 2 ) + A x p{x + y; ui,2u 2 ) 

+ C p(x; uji,2u 2 ) + C p(y; u)i,2w 2 ) + C 1 p{x]u 1 ,u] 2 ) + Cip(y; u u u 2 ), 
T(x, y) = 2(A p(x - y; 2ui,2u 2 ) + A p(x + y; 2u\,2u 2 ) 
+ A\p{x - y; ui, 2uj 2 ) + A x p{x + y;ui, 2u 2 )) 

x {Cop(x;ui,2uj 2 ) + C p(y; ui,2u 2 ) + C\p{x;u\, u 2 ) + C\p(y\ ui, lo 2 )) 
1 

- AAqCq p(x + ujj] 2a>i, 2uj 2 ) ■ p(y + ujj; 2wi, 2uj 2 ) 

j=0 

3 

- AA Ci ^2 P( x + Uj-, 2cji, 2w 2 )p(y + Wj-; 2wi, 2w 2 ) 

i=o 

- 4AiC p(x;u;i,2w 2 )p(y;u;i,2u; 2 ) 

1 

- AA1C1 ^2 P( x + w 2 fUJi,2oj 2 )p{y + uj 2j ;uJi,2uj 2 ). 

j=0 

(Trig-5 2 -S): ((sinh -2 Xt, sinh -2 2Af); (sinh -2 2At, sinh 2 2Ai}) 

-y(x) = Co sinh" 2 2Xx + Ci sinh 2 2Ax, w(y) = Cq sinh" 2 2Ay + C\ sinh 2 2Ay, 

u - (x — y) = Aq sinh -2 X(x — y) + A\ sinh" 2 2A(x — y), 

u + (x + y) = Aq sinh -2 X(x + y) + A\ sinh -2 2A(x + y), 

R(x, y) = Aq sinh -2 X(x + y) + Aq sinh -2 X(x — y) + A\ sinh -2 2A(x + y) 

+ Ai sinh -2 2A(x - y) + Cq sinh -2 2 Ax + Cq sinh -2 2Ay 

+ Ci sinh 2 2 Ax + C\ sinh 2 2 Ay, 
T(x, y) = 2{Aq sinh -2 A(x + y) + Aq sinh -2 A(x - y) + A\ sinh -2 2A(x + y) 

+ Ai sinh -2 2A(x - y)){C sinh -2 2Ax + Cq sinh -2 2 Ay + Ci sinh 2 2 Ax 

+ C\ sinh 2 2Ay) — j4oCo(sinh -2 Ax • sinh -2 Ay + cosh -2 Ax • cosh -2 Ay) 

- AAqCi (sinh 2 Ax + sinh 2 Ay + 2 sinh 2 Ax • sinh 2 Ay) 

- 4^iC sinh -2 2Ax • sinh -2 2Ay. 

(Rat-Ba-S): ((i -2 , t 2 ); (t~ 2 , t 2 )) 

v{x) = Cox' 2 + Cix 2 , w{y) = C y~ 2 + Ciy 2 , 

Aq Aq 

u~(x - y) = 7 + Ai(x- y) 2 , u + (x + y) = - — ■ — ^ + A\{x + y) 2 , 

(x-y) z {x + yy 



r^ + r^ + Mx-yf + Mx + yf + ^ + ^ + c^ + c.y 2 , 

(x — y) z (x + y) £ x z y z 



(x z - y 2 Y 
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(Toda-D^ } -S-bry): ((cosh2A, cosh4At); (sinh -2 2Xt), (sinh -2 2Xt)) 

v(x) = Co sinh -2 2Xx, w(y) = C\ sinh -2 2Ay, 

u~ [x — y) = Aq cosh2A(x — y) + A\ cosh4A(x — y), 

u + (x + y) = Aq cosh2A(x + y) + A± cosh4A(x + y), 

R(x, y) = Aq cosh 2A(x — y) + A cosh 2X(x + y) + A\ cosh 4A(x — y) 

+ Ai cosh 4A(x + y) + C sinh -2 2Xx + Ci sinh -2 2Ay, 
T(x,y) = 8Ai (C cosh 4Ay + Ci cosh 4Ax). 

(Toda-S^-S-bry): ((e -2Ai , e -4Ai ); (e 4At ), (sinh -2 2Xt)) 
v(x) = C e AXx , w(y) = d sinh -2 2Ay, 

u -(x -y) = A e~ 2X{x ~ y) + A ie - 4X(x - y \ u + (x + y) = Aoe~ 2X{x+y) + A ie - AX(x+y \ 
R{x,y) = A e- 2X(x - y) +A e- 2X{x+y) +A ie - 4X{x - y) +A ie - 4X{x+y) + 
T(x, y) = 4Ai (C cosh AXy + 2C x e~ AXx ) . 

(Toda-G^-S): ((e -2Ai ,e -4Ai ),0;(e 4Ai ),(e -4A <)) 
v(x)=C e AXx , w(y) = C ie - 4Xy , 

u~( x - y ) = Aoe~ 2X(x ~ y) + A ie - AX(x - y \ u + {x + y) = 0, 
R(x,y) = Aoe~ 2X{x - y) + A ie -4A(a:-J/) + C e AXx + C x e~ AXy , 
T(x,y) = 2A 1 (C e iXy + C ie - 4Xx ). 

(Trig-^-S-bry): ((sinh -2 At, sinh -2 2Xt), 0; (e -4At , e 4At )) 

v{x) = C e~ AXx + C ie iXx , w{y) = C e -4Ay + de 4 ^, 
u~~(x — y) = Aq sinh -2 X(x — y) + A\ sinh -2 2A(x — y), u + (x + y) = 0, 
R(x,y) = Aosmh' 2 X(x-y)+A 1 sinh -2 2A(x - y) + Coe~ AXx +C e~ AXy + de AXx +Cie AXy , 
T{x, y) = 2Ai sinh -2 2A(x - y){Coe~ AXx + Coe~ AXy + de AXx + de AXy ) 
+ 4A) sinh -2 A(x - y)(C e -2A(:!;+s/) + de 2 ^ x+y ^. 

(Rat d -£> 2 -S-bry): ((t 2 , t 4 ); (t -2 ), (t -2 )) 

= Cox -2 , = Ciy -2 , 

u~(x -y) = A {x - y) 2 + A 1 (x - y) A , u + (x + y) = A (x + yf + A\(x + y) A , 

R(x, y) = 2A (x 2 + y 2 ) + 2A x (x 4 + 6x 2 y 2 + y 4 ) + ^ + % 

T(x,y) = 32A 1 (C y 2 + C 1 x 2 ). 

Proof. (Ellip-i?2-S): We have the following from (|4.6j) . Lemma[4]and (12. 8j) . 

<2(p(i;wi,2u; 2 ); p(i;wi,2u; 2 )) = p(x; a>i, 2w 2 )p(y; wi, 2a> 2 ), 
<3(p(i; 2wi, 2ai 2 ); p(i; wi, 2u> 2 )) = Q(p(i; 2wi, 2w 2 ); p(i; 2u>i, 2w 2 ) 
+ p(t2Ai sinh -2 2A(x - y)(C + u\\ 2uj 1 , 2w 2 )) 

= p(x; 2wi, 2tv 2 )p(y; 2wi,2w 2 ) + p(x + 2a;i, 2w 2 )p(y + ur, 2u x , 2a> 2 ), 
Q(p(i; wi, 2w 2 ); p(t; ui, u 2 )) = p(x; uji,2co 2 )p(y, wi, 2w 2 ) 
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+ p(x +u 2 ;ui, 2u 2 )p(y + u 2 ; wi, 2w 2 ), 

3 

Q(p(t; 2ui, 2w 2 ); £>(*; vi,u 2 )) = ^ + uj; 2ui, 2u 2 )p(y + 2ui, 2u 2 ). 

3=0 

(Rat-B 2 -S) is given in [32^ (7.13)] but it is easy to check (|4.3p or prove the result as a limit 
of (Trig-BrS). (Rat d -.D 2 -S-bry) follows from T(t 2 ;t -2 ,0) = and T(t 4 ;t" 2 ,0) = 32y 2 (cf. Re- 
mark i)). Moreover (Trig-^-S), (Toda-D^-S-bry), (Toda-Cf } -S) and (Trig-Ai-S-bry) are 
obtained from the corresponding normal cases together with Lemma [U For example, Q for 
(Trig-£ 2 -S) is given by (f£7j) . and 

Q(sinh -2 At; sinh -2 2At) = Q ^sinh -2 At; i sinh -2 At - i cosh -2 At 

= - (sinh -2 At sinh -2 Ay + cosh -2 Ax cosh -2 Ay) , 

Q(sinh -2 2At; sinh -2 2At) = sinh -2 2Xx ■ sinh -2 2Ay, 
Q(sinh -2 2At; sinh 2 2At) = 0. 

Thus we get Theorem [3] from Theorem [2l 
4.3 Duality 

Definition 1 (Duality in B 2 , [23]). Under the coordinate transformation 

' x + y x — y s 



(x,y)^(X,Y) 



V2 ' V2 



the pair (P, P 2 — P 2 ) also satisfies (|4.ip . which we call the duality of the commuting differential 
operators of type B 2 . 

Denoting d x = d/dx, d y = d/dy and put 

L = P 2 ~P 2 - Q<9 2 -±tf + W - V Y + U~(d X + dyf + U + (d X - dy)'. 

Then the order of L is at most 2 and the second order term of L is 

-(«+ + u~ + v + w){d 2 + <9 2 ) - 2(u - - u+Jfl^ + 2wd 2 x + 2t;d 2 - (w - v){d 2 x - d 2 ) 

+ u -(d x + a y ) 2 + u+(d x - d y ) 2 = 0. 

Since L is self-adjoint, L is of order at most and the 0-th order term of L is 

--API + d 2 )(u + + u~ + v + w) + (u + + u - + u + wf - Avw - T - d x d y (u~ - u + ) 

d v)( w ~v) = (u + + u~ + v + w) 2 - Avw - T 
2 y 

and therefore we have the following proposition. 

Proposition 1 ( |23l I24j ). i) By the duality in Definition^ the pair (R(x,y),T(x,y)) changes 
into (R(x,y),T(x,y)) with 



R(x, y) = +w {^/f) + u+ (^ x ) + vr(y/2y), 
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ii) Combining the duality with the scaling map R(x,y) i— ► c~ 2 R(cx,cy), the following pair 
[R {x,y),T d {x,y)\ defines commuting differential operators if so does (R(x, y), T(x, y)j . This 
R d (x,y) is also called the dual of R(x,y), 

R d (x,y) = v(x + y) + - y) + ii + (2x) + u~(2y), 
T d (x, y) = R d (x, y) 2 - Av(x + y)w(x - y) - T(x + y,x -y). 
Remark 6. i) We list up the systems of type B2 given in Sections 14.11 and 14.21 

(u~(t),u + (t);v(t),w(t)) Symbo 

((p(i)); (p(t), pit + wi), p(t + w 2 ), p(t + w 3 )», (Ellip-B 2 

((p(i; 2wi,2u; 2 ), p(t; wi, 2w 2 )); (p(i;wi, 2w 2 ), p(t; w 2 )}), (Ellip-S 2 -S 

((sinh" 2 At); (sinh" 2 Xt, sinh" 2 2Ai, sinh 2 At, sinh 2 2Xt)), (Trig-B 2 

((sinh" 2 At, sinh" 2 2At); (sinh" 2 2At, sinh 2 2At)), (Trig-B 2 -S 

((t" 2 );(t" 2 ,t 2 , t 4 , t 6 )), (Rat-£ 2 

((t" 2 ,t 2 );(t" 2 ,t 2 )), (Rat-5 2 -S 

(cosh 2At); (sinh" 2 At, sinh" 2 2At), (sinh" 2 At, sinh" 2 2At)), (Toda-L^-bry 

(cosh At, cosh 2At) ; (sinh -2 At) , (sinh -2 At) ) , (Toda-D^-S-bry 

(e" 2A ');(e 2A *, e 4A '),(sinh" 2 At, sinh" 2 2At)), (Toda-fif } -bry 

(e" A ',e" 2A <); (e 2A <), (sinh" 2 At)), (Toda-B^-S-bry 

(e" 2Ai ),0; (e 2Ai ,e 4At ), (e" 2At , e" 4At )), (Toda-C^ 

((e- 2Ai ,e- 4Ai ),0; (e 4Ai ), (e" 4At )), (Toda-C^-S 

((sinh" 2 At), 0; (e" 2Ai , e" 4At , e 2Ai , e 4Ai )), (Trig-^-bry 

((sinh" 2 At, sinh" 2 2At),0; (e" 4At , e 4A *)), (Trig-Ai-S-bry 

((t" 2 ),0;(t,t 2 ,t 3 ,t 4 )), (Rat-^i-bry 

((t" 2 ),(t" 2 );(t 2 ,t 4 )). (Rat-D 2 -S-bry 

The first 6 cases above are classified by |24j as invariant systems. The systems such that at 
least two of u^, v and w are not entire are classified by [23]. (Trig-*) and (Toda-*) in the above 
are classified by [30] as certain systems with periodic potentials. 

We do not put ((t" 2 ,t 2 ),0; (t,t 2 )) in the list which corresponds to (Rat-A 1 -S-bry) because its 
dual defines a direct sum of trivial operators (cf. Section [9]). 

ii) Since 1 - sinh" 2 1 + 4 sinh" 2 2t = coth 2 t = t 2 - (2/3)t 4 + o(t 4 ) and sinh 2 t = t 2 + (2/3)t 4 + 
o(t 4 ), we have sinh 2 2Ax + sinh 2 2Ay - 2 coth 2 A(x - y) - 2 coth 2 X(x + y) = 8A 4 (x 2 + y 2 ) 2 + o(A 4 ). 
Hence the potential function 

%t + r Al ^ + ^ + ^ + C 1 (x 2 + y 2 ) +M x 2 + y 2 ) 2 



(x — y) 2 (x + y) 2 x 2 y 2 

is an analytic continuation of that of (Trig-£> 2 -S) but this is not a completely integrable potential 
function of type _B 2 

iii) The dual is indicated by superfix d . For example, the dual of (Ellip-B 2 ) is denoted by 
(Ellip rf -£> 2 ) whose potential function is 

3 

R(x, y) = Ap{2x) + Ap(2y) + C,(p{x -y + Uj) + p(x + y + uj)) 

3=0 
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and the dual of (Toda-Cj ) is 

«e- 2A ',e- 4A '), <e 2A V A *);0, (e" 4A *)) (Toda^-Cf ) 

since the dual of (u~(t), u + (t); v(t), w(t)^ is (w(t),v(t);u + (2t),u-(2t)). Similarly 

((p(i; u u 2uj 2 ), p(t; u u u 2 )); {p(2t; 2wi, 2^), p(2t; u x , 2cj 2 )}). (Ellip d -J3 2 -S) 

Since 4p(2t; 2a?i, 2W2) = p(i;^i,W2) etc, (Ellip rf -i?2-S0 coincides with (Ellip-B 2 -S) by replacing 
{ui,uj 2 ) by (2lo 2 , wi). 

Then we have the following diagrams and their duals, where the arrows with double lines 

represent specializations of parameters. For example, "Trig-i?2 => Trig-i?C2-reg" means that 
2 parameters (coupling constants C 2 and C3) out of 5 in the potential function (Trig-^) are 
specialized to get the potential function (Trig-iJCVreg) with 3 parameters. For the normal cases 
see Definition [5] and the diagrams in the last part of the next section (type B n ). 

Hierarchy of Normal Integrable Potentials of type B 2 

Ti\g-BC 2 -xeg — > Toda-Z?2-bry 

C 2 =C 3 =0 ft 5:3 C =C 1= ft 5:3 

Trig-^2 Toda-S^-bry U Toda-5^ 

/ / I 

Ellip-5 2 -> Toda-D^-bry Toda-C^ ™ Toda-BC 2 

c =c 1= o 

\ 

5 - 3 

Trig-i?2 Trig-^i-bry => Trig-Ai-bry-reg 

c 2 =c 3 =o 

\ \ 

Rat-i? 2 —* Rat-^i-bry 



Hierarchy of Special Integrable Potentials of type B 2 

TrigW-B 2 -S-reg -> TodaW-D 2 -S-bry Rat d -D 2 -S-bry 

Ci=0 ft 4:3 C =0 ft 4:3 / 

TrigW-5 2 -S -► Toda^-B^-S-bry ^ Tbda^-B^-S 

/ / I 

Ellip-B 2 -S -> TodaW-D^-S-bry TodaW-C^-S ^ Toda^-^-S 

Co=0 

\ 

Trig( d )-B 2 -S -> Trig( d )-^i-S-bry ^ TrigW-Ai-S-bry-reg 

\ 

Rat-5 2 -S 

Definition 2. We define some potential functions as specializations. 

(Trig-i?2-S-reg) Trigonometric special potential of type B 2 with regular boundary conditions is 
(Trig-Sa-S) with d = 0. 

(Toda-Z?2-S-bry) Toda special potential of type D 2 with boundary conditions is (Toda-B^-S-bry) 
with Co = 0. 

(Toda-i^l -S) Toda special potential of type B^ is (Toda-U^-S-bry) with C\ = 0. 
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(Toda-i? 2 -S) Toda special potential of type -B 2 is (Toda-C^-S-bry) with Cq = 0. 

(Trig-^4i-S-bry-reg) Trigonometric special potential of type A\ with regular boundary conditions 
is (Trig-^-S-bry) with d = 0. 

Remark 7. We have some equivalences as follows: 

(Ellip-£ 2 -S) = (Ellip d -£ 2 -S), 
(Rat-5 2 -S) = (Rat d -B 2 -S), 
(Trig-SCs-reg) = (Trig^-S-reg), 
(Trig-^i-bry-reg) = (Toda d -Z) 2 -S-bry), 
(Toda-D 2 -bry) = (Trig d -,4i-S-bry-reg), 
(Toda-B^) = (Toda d -^ 1} -S), 
(Toda-SCa) = (Toda d -£ 2 -S). 

5 Type B n (n > 3) 

In this section we construct integrals of the completely integrable systems of type B n appearing 
in the following diagram. The diagram is also given in P Figure ELI], where (Toda-S^-bry) 
is missing. The most general system (Ellip-i? n ) is called Inozemtzev model (cf. [12]). 

Hierarchy of Integrable Potentials with 5 parameters (n > 2) 

Toda-D^-bry -> Toda-i^-bry -» Toda-C^ 

Ellip-B n -> Trig-S n -> R&t-B n 

\ \ 

Trig-^n.i-bry ->• Rat-^ n _i-bry 

5.1 Integrable potentials 

Definition 3. The potential functions R(x) of (II. ip are as follows: 
Here A, Co, C\, C 2 and C3 are any complex numbers. 

(Ellip-B n ) Elliptic potential of type B n : 

n 3 

A(p{xi - Xj-;2wi,2w 2 ) + p(xi + Xj-;2o;i,2a; 2 )) + y^^Cjp(x k + Uj\ 2a>i, 2w 2 ), 

l<j<j<n fc=l j=0 

(Trig-i? n ) Trigonometric potential of type B n : 

j4(sinh -2 A(xj — Xj) + sinh -2 A(xj + Xj)) 

l<i<j<n 

E/ 2 2 2 C3 2 1 

I Co sinh \xk + Ci cosh Ax^ + C 2 sinh Xxk + — sinh 2Ax/c I , 

k=l ^ ' 
(Rat--B n ) Rational potential of type B n : 

E ( (x .- x .)2 + (x . t x v ) + 2 + c ^ + + 
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(Trig-A n _i-bry) Trigonometric potential of type A n _\ with boundary conditions: 

n 

Asinh" 2 X( Xi - X j) + J2{C e~ 2Xxk + de~ 4Xxk + C 2 e 2Xxk + C 3 e 4Xxk ), 

l<i<j<n k=l 

(Toda-S^^-bry) Toda potential of type Bn^ with boundary conditions: 
n-l 

^ Ae -2\(x t -x t+l ) + Ae -2\{x n _^ n ) + C Qe 2\x 1 +C7ie 4Ax 1 ^ ^-2 A^ + Cg skill" 2 2Ax„, 
i=l 

(Toda-C^) Toda potential of type C^ 1 '- 

71-1 

J2 Ae~ 2X(xi - Xi+l) + C e 2Xxi + C ie 4Xxi + C 2 e~ 2Xxn + C 3 e~ 4Xxn , 
i=i 

(Toda-D^-hry) Toda potential of type Dn^ with boundary conditions: 

n-l 

Ae -2\(xi-x i+1 ) _|_ Ae -2X(x n -!+x n ) _|_ Ae 2\(x 1 +x 2 ) 



i=l 



+ Cosinh 2 Axi + C\ sinh 2 2Axi + C2 sum 2 Ax n + C3 sinh 2 2Ax n , 



(Rat-^4 n _i-bry) Rational potential of type A n ^\ with boundary conditions: 



A n 

Yl (x . _ x .\2 + I]( g ox fe + Ci4 + c 2 x\ + c 3 4). 

l<i<j<n i> k=l 



Remark 8. In these cases the Schrodinger operator P is of the form 

1 n rft 



j=l 3 



3 n 
R{x)= Yl ( U *i-ej( X ) + U e l +e j {x)) +^2C j V 3 (x), V ] (s) = ^ ^ fc (x) . 

i<i<j<n i=o fc=i 

Here 

<9 a n Q (x) = d^v^x) = if a, i G 1" satisfy (a, a) = (b,f3) = 0. 

The complete integrability of the invariant systems (Ellip-fi„), (Trig-S n ) and (R&t-B n ) is 
established by [29] . We review their integrals and then we prove that the other five systems are 
also completely integrable by constructing enough integrals, which is announced by [30J. The 
complete integrability of (Trig-^n^i-bry), (Toda-C,^), (Toda-D^-hry) and (Rat-^4 n _i-bry) is 
presented as an unknown problem by [8] and then that of (Toda-S^^-bry), (Toda-C^) and 
(Toda-D^-bry) are established by [TBI 021 EJJ using i?-matrix method. The compete integra- 
bility of (Trig-j4 n _i-bry) and (Rat-y4 n _i-bry) seems to have not been proved. 

Definition 4 ([25} 129]). Let u a (x) and T°(v J ) are functions given for a S £(-D n ) and subsets 
I = . . . , ik} C {1, . . . ,n} such that 

Uaix) = and d y u a = for y G W n with {ct,y) = 0. 
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Define a differential operator 
n 1 

P n{C) = k w n _ k )\ E {l{w(l),-Mk)}(. C ) • A \w(k+l),...,w(n)}) 

k=Q ' V '* UlGS n 

by 

&{h,...,i h } = E 2 k j\(k-2j)l (5,1 ^ 



o<i<[|] 
x 



E £ ( W ) W ( U ^ 1 ~e l2 (x)u ei ,^ eH (x) • • • % 3 r e, 2 . (x) ■ ^y+i^fey+j, • ■ ■ 

^ 1 ,..,u-}( C ) = E E (5.2) 

„=i / 1 u--n/„={u,...,i fe } 

Tj = {-I)* 1 ' 1 (est - J2 CjTfW)) > (5-3) 

where 

S {h,i 2 ,...,i k } = 2 E Uj K n - e!2 (^)« ei2 -e l 3(^---u eife _ 1 - eifc (^)), 
i«6lV(B fe ) 

S° = 0, = 1, S° {IJ} = 2u ei _ ej (x) + 2n ei+ej (x), 
W) = 0, T { ° fc} (^) = 2<(x) for l<A:<n, 
Q0 = 1, = g {il j 2} = T {il} T { , h} + T {ili , /2} , 

In the above, we identify W(-Bfc) with the reflection group generated by w eik and w e ^_ ei ^ i 
(z/ = 1, . . . , k — 1). The sum in (|5,2p runs over all the partitions of the set / and the order of the 
subsets Ii,...,I u is ignored. 

Replacing dj by £j for j = l,...,n in the definition of Afj,^^! and P n (C), we define 
functions A/^ ^ i and P n (C) of (x,£), respectively. 

We will define u a {x) and Tf{v ] ) so that 

[P n (C),P n (C')] =0 for C,C'eC. (5.4) 
Then putting 



q°i = qi 



C*=0' 



P n - P n (0) - ^ fc w _ j^j E (^(l),..,w(fc)} A {w(fe+l),..,w(n)})> ( 5 - 5 ) 

fc=o ' v ; ' wee n 

^"EEjia^iLMi E E ••• 5 '/iC({i+i r ..,fe}) A H{fe+i,...,n})' ( 5 - 6 ) 

i=j k=i V Jy ' u)€6„ /iU---n/j=to({l,...,i}) 

we have P n (C) = ]T C J 'P n -j and CCU and then 

[P,P,]=0 for l<i<j<n, 

p i n q2 3 

-Y = -2Ep + E K i -e J (x)+U ei+e ,(x))+^jC J V(x). (5.7) 
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Remark 9. i) When n = 2, T(x,y) in the last section corresponds to T± 2 , namely 
T(xi,x 2 ) = T 12 \ c=(j . 

ii) Put 

n 

U{x)= ^ ( U ij( X t - X j) + U tj( X i + X j)) and V ( X ) = ^2v k (x k ) 

l<i<j<n k=l 

in (|1.3j) and let Tj(U ; V) be the corresponding Tj given by (|5.3p . Then [29l Remark 4.3] says 
T/(co^;ciy + c 2 W) = c* I - 1 c 1 T I (U;V)+c* I - 1 c 2 T I (U;W) for cj G C. (5.8) 

iii) The definition (|5.3|) may be replaced by 
#i 



/ #1 3 \ 

T 7 = (-I)*'" 1 £ (-2A)^ 1 (z, - 1)! • 5 7 V • • SI - CjTj> j ) 

V is=i hn - ui„=i j=o / 



(5.9) 



because A can be any complex number in |29} Lemma 5.2 ii)] when v = C. Note that we fixed 
A = 1 in [29] . Combining (|5,9p and (|5.8p , we may put 

tj = (-if 1 - 1 (cs? + C 'Y^ C- 1 y, s °h • • • 5 t " E C i T ° (W) 

V u>2 hU~llI v =I j=o J 

for any c, d E C and hence 



Tj = (-I)* 1 ' 1 CS] + j> E s l - E (■ 

V i/>2 / 1 u--n/„=7 j=o / 

for any C2, C3, . . . G C. 

Theorem 4 (Ellip-S n , [25], [291 Theorem 7.2]). Put 

Ue i ±e j (x) = Ap (xi ±Xj;2ux,2u 2 ) for 1 < i < j < n, 

v 3 ek (x) = po(xk + u)j] 2u\, 2to 2 ) for 1 < k < n and < j < 3 

and 

#1 

W) = E E ("^"V - 1)! • S h (v 3 ) ■ ■ ■ (5.10) 

S{ tl ,...,i k }{v J ) = Yl ^(e il )( X )^(en-e« 2 )( a; )^(e i2 -e i 3)(^)--^(e ifc _ 1 - eife )(^)- (5-11) 
T/ien (ED feo/ds. 



3 

Example 1. Put v 3 k = vi k ,Vk= J2 Cj v l' and w ij = u ij ^ u tj- Then 

A {i,2} = did 2 + iij~ 2 - = did 2 + w7/ 2 , 
A {i,2,3} = <9i<9 2 <93 + u>7/ 2 <9 3 + w^dx + 
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A{i,2,3,4} = did 2 d'id A + w- iA did 2 + w 24 did 3 + w 22t did A + w u d 2 d 3 + w 13 d 2 d 4 
S 1 ^)=2v{, 

<%,2}(^) = 2vi U l2 + 2v l U 12 + 2v 2 U 12 + 2v 2 U 12 = 2 ( V 1 + v i)( U 12 + "12)1 
^{l.^S}^) = 2l, l u r2 u 23 + 2v { u 12 U 22, + 2v \^X2 U 2Z + 2v L u 12 U 23 + ' ' ' 

= 2{v{ + v^)wf 2 w^ 3 + 2{v{ + vQw&wfa + 2(v J 2 + v^)wf 2 wf 3 , 
3 

r {1} = - £ c/r { ° 1} (^) = c - 2%, 

i=o 

3 

3=0 

^{1} = T {1}' 

9{1,2} = T{l}T{2} + T{1,2}-, 

9{i,2,3} = T{i} T {2}T{z} + r{ li2 }T{ 3 } + T{ li3 }T{ 2 } + T{i}T{ 2j 3} + 2{ 1)2i3 }. 
If T°(v*) and S^,...,^^') are given by ([530]) and ([5TTTT) . then 
T{ 1} (vi) = 2v{, 

T^ 2} (v j ) = S {1 #(v>) - AT° {1} T° {2} = 2(v{ + «>+ - AAv{v{, 

r { i, 2> 3}(^) = S{i,m&) ~ 2 M4S{2M vj ) + 4s { iM vj ) + ^{i,2}(^')) + i6A 2 444. 

In particular, if n = 2, then 

P 2 (C) = Aj lj3} + <?{i}Aj 2} + |?{2}A{ 1} + g { i i2 } = {did 2 + u{ 2 - u\ 2 ) 2 

3 

+ T {1] d 2 2 + T {2} d\ + T {1} T {2} - CS° {h2} + C 3 T° {1 ^) 

3=0 

= (did 2 + u\ 2 - u\ 2 f + {C- 2vx)dj + (C - 2v 2 )d\ + (C - 2^)(C - 26 2 ) 

3 

- 2C(uf 2 + u+ ) + 2(ui + #2)^2 + u+ 2 ) - AA ^ C jV {v j 2 = C 2 - 2P ■ C + P 2 



3=0 



with 



P = ~\( d l + d i) +vi+v 2 + + u+ , 

3 

P2 = {d^+u^-uf^ 2 -2vid 2 -2v 2 dl+4v 1 v 2 + 2(v 1 +v 2 ){u^ 2 +uf 2 )-4A^2c j viv J 2 , 

3=0 

which should be compared with (|4.2p . (|4.4p and (|4.6p . 
In general 

n n 

* = £(*?*}+«?*>)- £ sfo = - ^ - 2 £ 

fc=l l<i<j <n. k=l l<i<j'<n 

ft- E E i>t.}* 2 w+ E 

Ki<j<n l<i<n j=l l<«<3<n 

l<j<n,v£j 
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E S {i,j}( A {k} + 1°{k}) + E S {i,j} S {k,£} + E ^{ij.fe} 

l<i<j<fc<n 



l<i<j'<n 
l<fc<n. kj^i, j 



l<i<j<n 
i<k<£<n 



= + w~ 3 ) 2 - E ^9] + E 45^ + ^ C.Tj, 

-2£<4 - 2%) + 4 £ iu±u;+ + 2 E «,+«,+ . 
Here if f|5. 10|) and (|5. 1 1 j) are valid, then 

T° {kA {v 3 ) = 2(4 + 4)w+ - 4Av{vi 
The commuting operator P3 of the 6-th order is 

P 3 = E A hM + E «h A h} + E «?U*> - E S lkA* 

- E 5 {M} ? » A {j} ~ E S {k,i} q kj} + E 5 '{i > i,fe,Am} + E S '{i,i,*:/} 5 '{/x^} 
1 \ A c*o c*o I \ A no no c*o j_ \ A qo qo qo qo 

In Theorem H] the Schrodinger operator is 



(5.12) 



(5.13) 



,2 



P = --E^T + A E (P( x * - x i) + P( x i + + E C J E + 

fe=0 Xk l<i<j<n j=0 k=l 

and the operator P2 satisfying [P, P2] = is given by (|5. 12 j) and (|5.13p with 

3 

Vk = ^2c„p(x k +u v ), v{ = p(x k wf j =A(p(x i -Xj)±p(x i + Xj)). 

5.2 Analytic continuation of integrals 
Theorem 5 (Toda-.DW-bry). For 



U ei+e] [x) = \ 



Ae -2X( Xl -x i+1 ) (j = i+ 1) ; 

(|i-i|>l), 
' Ae 2X(x 1+ x 2 ) (i+j = 3), 

Ae -2\(x„- 1+ x n ) (j + j = 2n - 1), 
.0 (i + j i {3,2n-l}), 

v\{x) = 5\k sinh" 2 \x\, v k( x ) = °~ik sinh" 2 2Axi, 
v 2 k (x) = 5 nk sinh" 2 Xx n , v\(x) = 5 nk sinh" 2 2Xx n , 

we have commuting integrals Pj by (|5.5p . ()5.ip . (|5.2p . (|5.3p and 



(5.14) 



S 



{fc} 



1 for 1 < k < n, 



Sf = if I ^{k,k + !,...,£} for l<k<£<n, 



S 



{k,k+l,...,e} 



2A £-k+l( e -2\(x k -xt) + e 2A(x 1+ x,) + Sine -2X(x k +x n )^ 



Tf k} (vi) = 2v{(x) for 0<j<3, fc = l,...,n, 
lf(v°) = if I^{l,...,k} for k = l,...,n, 
T 7 V)=0 if I^{k,...,n} for k = l,...,n, 
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T,V) = 7/V) = if #/ > 1, 

T {i,...,fc}( w °) = %A k ~\e 2Xx * + 5 kn e- 2Xx ") for k > 2, 

T? n - k+1 ,..., n} (v 2 ) = SA^ie-^-k+x + S kn e 2Xxi ) for k > 2. 



Proof. Put 



1-1 k- 1 n-1 

Xi -U> 2 , • • • , Xk 7^2, ■ ■ ■ ,X n -U) 2 

n — 1 n — 1 n — 1 

e 2Aw2/(n-l) / j _ 1 / 7 — 1 

"e^e,- (5) = A ^ P0 ( Xi - - - ^ U) 2 T i x j — - - - UJ 2 ) \ 2(jJ±, 2uJ 2 



^fc(^) = — T9 - Po ( Xk -u>2 + Wj-; 2u?i, 2a?2 for < j < 3, 1 < fc < n. 



(-1)' / fe — 1 
W = Pol a* - —[^ 

When uj 2 — > oo, u eiTej (x) and «| (i = 0,1,2,3) converge to u eiTej (x) in (|5.14p and 

v k (x) = 5ik sinh -2 \x\, v k {x) = 8\ k cosh -2 Axi, 
vl(x) = 5 nk sinh" 2 Ax n , v k (x) = 5 nk cosh -2 Ax n , 

respectively. Under the notation in Theorem HI let Si(v e ) and T°(v ) be the functions defined 
in the same way as Si(v ) and Tf(v e ), respectively, where (u eiTej (x), v e k (x)) are replaced by 
(u ei zp ej (x), v k {x)). Then by taking the limits for uj 2 — > oo, Tfffi) converge to the following Tf(v e ) 

ff(v°) = ff(v l ) = if I^{l,...,k} for k = l,...,n, 
ff(v 2 ) = TjV) = if I^{k,...,n} for k = 1, . . . , n. 

If k > 2, then 

# 7 f „2Aw 2 /(n-l) \ 

v=\ i 1 n-m„=i 1 



#£ / 2Aw 2 /(n-l) \ 

,/— 1 r. tt . it r — r V / 



'„~.0\ i- S . /-On ^ e 0, /~0> 



= I™ 5 { i ,...,&}(£ ) - lim 5m (v )A — ^ «S , / 2) ... ) a(«", 

= 2A fc " 1 sinh- 2 Ax 1 (e- 2A ( a;i - :Efc) + e 2A ( Xl+a;fc ) + fl^e 2 ^* 1- *") + £ fcn e~ 2A(a:i+a:Tl) ) 
- 2 sinh" 2 Ax! • 2A k ~ 1 {e 2Xx * + <5 fcn e- 2Ax ") = SA*- 1 ^ 2 *** + b kn e~ 2Xxn ), 
T{i,...,k}( v ^ = 2Ak ~ X cosh- 2 Ax 1 (e- 2A(a;i - Xfc) + e 2A(:ri+Xfc) + ^e^^ 1 -^^ Jfcne"^^ 1 ^")) 

+ \A k ~ x cosh" 2 Axi(e 2A:Cfc + 5 kn e- 2Xxn ) = %A k - x (e 2Xxk + 8 kn e~ 2Xxn ), 
T° {n . k+1 ,..., n} (v 2 ) = 8A k -\e- 2Xx -^+5 kn e 2Xx ^, 
T° {n . k+1 ,..., n} (v 3 ) = 8 A k -\e- 2Xx -^ + 5 kn e 2Xx -). 



Replacing v 1 and v 3 by (l/4)(t>° — v 1 ) and (l/4)(v 2 — u 3 ), respectively, we have the theorem by 
the analytic continuation given in Lemma [2j ■ 

As is proved by [29], suitable limits of the functions in Theorem [J] give the following theorem. 

Theorem 6 (Trig-S„, [29J Proposition 6.1]). For complex numbers X, C , Co, . . ., C3 and A 
with A 7^ 0, we have (|5.7p by putting 

Uei±eAx) = A sinh" 2 A(xj ± x k ), 
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v e k ( x ) = sinli 2 Xxk, v l k ( x ) = cosh. 2 \Xk, 
v e k ( x ) = sinh 2 Xxk-, v e k ( x ) = 7 sinh 2 2Xx k 

and 

Tj = (-lf^iCS ! - C T?(v°) - dTfiv 1 ) - CzS^v 2 ) - G 3 Sj(v 3 ) 

+ 2C 3 (S^Sj^ + Sj^S^ + SISj^v 2 ))), 

hui 2 =i 
#i 

T ?(v°) = E E (-^rV - • • • s/>°), 

r /V) = £ £ - 1)!*!^ 1 ) • • • Si„(v l ). 

Theorem 7 (Trig-A^-i-bry). For 

u ei - e] (x) = A sinh -2 A(x; - xj), u ei+e .(x) = 0, 

< (x) = e" 2A ^ , < (x) = e" 4A ^ , < (x) = e 2A ^ , < (x) = e 4 



we /ioue (|5.7|) by putting 



Tj = (-l f 1 - 1 (cS°j -^C^ 7 (<) + 2^T (^S^V^V^S/^ 2 )^ 2 ))) . 
Proof. Putting 

i±e . = ^sinh" 2 A((xi + JV) ± (ay + TV)), {}° = - e 2XN sinh" 2 A(x fc + N), 



e,_i_cj 



ul = i e 4A7V sinh" 2 2A(x fc + JV) = — e 4A7V (sinlr 2 2A(x fc + AO - cosh" 2 2A(x fc + AT)). 
4 16 

v 2 = Ae~ 2XN sinh 2 A(x fc + N), v 3 k = 4e~ 4XN sinh 2 2A(x fc + N), 
x = (xi + N, x 2 + N, . . . , x n + N) 

under the notation in Theorem [6l we have 

(u ei ^ e u ei+e vl,vl,vl,vl) := lim (u Ci _ e ,u ei+e vl,vl,vl,vl) 

N— too 

= (Asmh~ 2 X( Xi - x k ),0,e- 2Xxk ,e~ 4Xxk e 2Xxk ,e AXxk ), 
lim i e 2AAr T / ( W °)(x) = 5 / (^ ), 



l im l e ^(2y(«»)(£)-2y(i; 1 )(£)) = 5r(t; 1 )-2 £ S h (v°)Sj 2 (i 



lim 4e- 2A7V T / °(^ 2 )(x) = 5/(v 2 ), 

7V-»oo 

lim 4e- 2AAr T/(t; 3 )(x) = 5 / ^ 3 )-2 £ S/^ 2 )^ 2 ). 

hlll 2 =l 

Here Sj(v e ) are defined by (|5.11|) with u ei ± ej and tc, replaced by u ei ± e . and tig , respectively. 
Then the theorem is clear. ■ 
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Theorem 8 (Toda-fi^'-bry). For the potential function defined by 



Uei+eAx) 



' Ae- 2X{ - Xi ~ Xi +^ if j = i + l, 
if l<i<i+l<j<n, 

^ e -2A(x n _ 1 +a; n ) -j i = n - 1, j = n, 
if 1 < i < j < n and i / n — 1 , 

v° k {x)=5 kl e 2Xx \ vl(x) = 5 kl e AXx \ 
v 2 k {x) = 5 kn sinh" 2 Xx n , v\(x) = 5 kn sinh" 2 2\x n , 

we have (|5.7|) with 

S° k} = 1 for l<k<n, 

S? = if I ^ {k,k + !,...,£} for 1 < k < t < n, 



S {k,k+i,...,£} = 2A e - k+1 {e- 2X ^- x ^ + 5 £n e- 2X ^ +x ^), 
lf k} (v j ) = 2vi(x) for 0<j<3, k = l,...,n, 
lf(v°) = if I^{l,...,k} for k = l,...,n, 
T?(v 2 ) = if I^{k,...,n} for k = l,...,n, 
T?(v 1 )=T?(v 3 ) = if #/>l, 



T^... M (v U ) = 2A k -'(e 2 ^+5 kn e 



T 



{U+l,...,n}^ ) = &4 



fc-1 -2Ax 



Proof. Suppose Re A > 0. In (Toda-D^-hry) put 

x, = (x± — {n — 1)N, . . . , x k — (n — k)N, ... ,x n — (n — n)N), 

' A e -2\N e -2\{xi-{n-i)N-x l+l + {n-i-l)N) f j _ ^ _|_ ^\ 

.0 (|j-i|>l), 
'^4 e -2AAr e 2A(a:i-(n-l)Af+X2-(n-2)Af) ^ + j _ 3^ 
^ e -2A7V e -2A(x n _ 1 -JV+x n ) (% + j = 2n - 1), 

.0 (* + i^{3,2n-l}), 

e 2A(n-l)jV 



U, 



Ud+ej 



V k = Sl k - 



sinh"^A(rEi - (n - 1)N), 



e 4\(n-l)N 

4 = $lk — t 



sinh"^2A(xi - (n - 1)N), 



v k = o~nk sinh Xx n , v k = 5 nk sinh 2Xx r 



(5.15) 



and we have (|5.15p by the limit N — > oo. Moreover for k > 2, it follows from Theorem [5] 
and (ED that 



T {1 ,..., k} (v 1 )=f {1 _ k} (i 3 )=0, 

f {1 ,..., k} (v°) = I e 2A(n-l)7V (Ae -2A7V )fc -l (8e 2A(x fc -(n- fc )7V) + g^-aAz^ 



{n-/c+l, 



= 2A fc " 1 (e 2Axfc +4„e" 2Ax "), 

>n} (i) 2 ) = ( J 4 e -2AJV ) fc-l ( g e -2A( !En _ fc+1 -(fe-l)iV) +8(5ifce 2A(x 1 -(n~l)JV) ) 



|fc-l„-2Ax 



n-k+l 



which implies the theorem. 
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Theorem 9 (Toda-C^ 1 )). For the potential function defined by 



U ei -ej [X 



f Ae -2X(xi-x i+1 ) if j = i + l, 

[0 if \<i<i + l<j<n, 

u ei+e .(x) = for 1 < i < j < n, 
v° k (x) = 5 kl e 2Xx \ v\{x) = 5 kl e iXx \ 
v 2 k (x) = 5 kn e~ 2Xx % v 3 k (x) = 5 kn e~ AXx \ 

we have (15.71) with 



S° {k} = 1 for l<k<n, 

Sf = if I ^{k,k + l,...,£} for l<k<£<n, 

no _ 9 \l~k+l ~-2X(x k -xg) 

D {k,k+1,...,£} — zyl e > 

T?(v°) = i 
T?(v 2 ) = i 
T?(v 1 )=T?(v 3 ) 



T 



for 


< j < 3, jfe 


= l,...,n, 




. . . , fc} /or 


k = 1, . . . ,n. 


I^{k, 


. . . , n} /or 


k = 1, . . . , n 


if 


#/> 1, 




-l e 2Xx k 


/or fe > 2, 




2A k - 1 e- 2Xx "- k +i for 


k>2. 



\n-k+l,...,n}( V ) — 

Proof. Substituting x k by x k + R for k = 1, . . . , n and multiplying i>P, ui, i> 2 and i>| by e~ 2A ^, 
g -4A_R^ (1/4^2^^ anc [ (l/4)e 4Afi , respectively, we have the claim from Theorem[8l ■ 



Ue z -e 3 {X) = — — 2 , We l+ e,(x) = 0, v£(s) = 4 + \ 



Theorem 10 (Rat-A„_i-bry). We We (J5T7D if 

A 

(Xi Xj)' 

T I = (-l)* I ~ 1 (cS° I -^2c j S I ^ j )+ E ^(S^^ + S?^ (*/>)) 

V j=0 /iU/ 2 =-f 

+ E G3(5/ 1 ( W 1 )S / o 2 + S /l ( w )^ 2 ( v ) + S / o 1 5 /2 (^ 1 ))). 

/iU7 2 =/ / 

Proof. Put 

n ei _ ej = A 2 sinh" 2 A(xj - Xj), u^+e^ = 0, 

vl = ^(e 2Xx * - 1), v\ = ^(e 2Xx * + e~ 2Xx - - 2), 

£ 2 = J-( e ^k_ 3e 2Xx k _ e -2Xx k ~3 = _l_( e 4Ax fe+e -4Ax fc _ 4e 2Ax fe _ 4e -2A :Cfc+6 ^ 

Then taking A -> we have the required potential function. 
Owing to (Trig-A n _i-bry) and Remark we have 



A->0 

1 2 



A ^ (E g2AXfc ) (E e2AXfc ) - 4 ^ft° 2 
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= l(s h (J2 *k) S° h + S° h S l2 (J2 *k) ) , 

ft a2 iw {s h (E e 2XXk ) ft* (E e2A i + fti (E e ~ 2A i ^(E e_2A i - 8 ^A 



^ (X (E **) ^ + ^ (E **) (E x >) + ^ (E 



Xk 



and thus the theorem. I 
As is proved by [22] , suitable limits of the functions in Theorem [U] give the following theorem. 
Theorem 11 (Rat-S„, [291 Proposition 6.3]). Put 
A A 



u ei-eA x ) — , \ 2 ' u ej+ej 0*0 



(x^ Xj^ {xi -\- Xj^ 

v K x ) = x k 2 > v U x ) = x h v U x ) = x t v li x ) = x t- 

Then (|5.7p holds with 

Tj = (cS 7 ° - C T?(v°) - dS^v 1 ) - CzS^v 2 ) 

+ 2C 2 E (S h (v 1 )S] 2 + SlS h (v 1 ))-2C 3 SKv 3 ) 

hUI 2 =I 

+ C 3 E ("5/ 1 (^ 1 )S /2 (^ 1 ) + 25' Jl (« 2 )5? a + 25y i 5 /2 (u 2 )) 
/in/ 2 =/ 

- 24C 3 E ^(v 1 )^, + Sf^iv 1 )^ + SlSlS^v 1 ))] , 
/1 11/211/3=/ / 

#1 

T ?(v°) = E E (-^)"~V - !) ! • • • • sua 

i/=i /!u-n/„=/ 

Definition 5. We define some potential functions as specializations of potential functions in 
Definition [3j 

(Trig-A„_i-bry-reg) Trigonometric potential of type A„_i with regular boundary conditions is 
(Trig-^ n _!-bry) with C 2 = C 3 = 0. 

(Trig-A n _i) Trigonometric potential of type A n -\ is (Trig-A n _x-bry) with Co = C\ = C 2 = 
C 3 = 0. 

(Trig-.BC n -reg) Trigonometric potential of type BC n with regular boundary conditions is (Trig- 
B n ) with C 2 = C 3 = 0. 

(Toda-D n -bry) Toda potential of type D n with boundary conditions is (Toda-B^ -bry) with 
C = Ci = 0. 

(Toda- i?! 1 '* ) Toda potential of type Bn^ is (Toda-Si^-hry) with C 2 = C3 = 0. 
(Toda-Z?i 1 ' ) ) Toda potential of type is (Toda-Z^-hry) with Co = C\ = C 2 = C3 = 0. 
(Toda-y4 n _i) Toda potential of type A n _\ is (Toda-C^) with Co = C\ = C 2 = C3 = 0. 
(Toda- BC n ) Toda potential of type B n is (Toda-Cjj 1 ^) with Co = C\ = 0. 
(Ellip-D n ) Elliptic potential of type D n is (Ellip-i? n ) with Co = C\ = C 2 = C3 = 0. 
(Trig-D n ) Trigonometric potential of type D n is (Trig-£>„) with Co = C\ = C 2 = C3 = 0. 
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(Rat--D n ) Rational potential of type D n is (Rat--B n ) with Co = C\ = C2 = C3 = 0. 

(Toda-D n ) Toda potential of type D n is (Toda--Bn^-bry) with Co = C\ = C2 = C3 = 0. 

(R&t-B n -2) Rational potential of type B n -2 is (Rat-f? n ) with C2 = C3 = 0. 

(Rat-j4 n _i-bry2) Rational potential of type A n _\ with 2-boundary conditions is (Rat-A„_i-bry) 
with C2 = C3 = 0. In this case, we may assume Co = or C\ = by the transformation 
x k l— > x k + c (k = 1, . . . , n) with a suitable c G C. 

Then we have the following diagrams for n > 3. Note that we don't write all the arrows in 

5'3 

the diagrams (ex. (Toda-Z? n -bry) — > (Tod&-BC n )) and the meaning of the symbol is same 
as in the diagram for type B2. Namely, 5 parameters (coupling constant) in the potential are 
reduced to 3 parameters by a certain restriction. 



Hierarchy of Elliptic- Trigonometric- Rational Integrable Potentials 

Rat-J3 n -2 Ellip-A, 

ft5:3 I 

Rat-£ n \ Trig-£> n -> Rat-D„ 

T U.i 

Ellip-£ n Trig-B n ^ Trig-SC n -reg Ellip-A^ 

Trig--4n-i-bry ^$ Trig-^ n _!-bry-reg %> Trig-^ n _i 

5'3 3'1 

Rat-^4 n _i-bry =4> Rat-A n _i-bry2 ^ Rat-A„_i 



Hierarchy of Toda Integrable Potentials 

31 

Trig-f?C ra -reg — > Toda-D n -bry Toda-D n 

fl~5:3 1>5:3 1>3:1 

Trig-5„ -> Toda-^ 1} -bry Toda-5^ 

/ /* \ 

Ellip-5,, Toda-D^-bry U Toda-D^ Toda-C^ 

^5:1 / ^5:3 

Ellip-L> n Trig-^4 n _i Toda-5C n 

/ \ ^3:1 

Elhp-A re _i -> Toda-A^ -> Toda-A n _i 



6 Type £> n (n > 3) 

Theorem 12 (Type -D n ). T/ie Schrddinger operators (EHip-D n ), (Trig-D n ), (Rat-D n ), (Toda- 
D ( n ] ) and (Toda-L> n ) are in the commutative algebra of differential operators generated by 
Pi, P2, . . . , P n -\ and A{ l5 n \ which are the corresponding operators for (Ellip-i? n ), (Trig-i? n ), 

(Rat-B rt ), (Toda-i^-hry), (Toda-D n -bry) with Co = C\ = C2 = C3 = 0, respectively. 

Proof. This theorem is proved by [29] in the cases (Ellip-D n ), (Trig-D n ), (Rat-D n ). Other two 
cases have been defined by suitable analytic continuation and therefore the claim is clear. ■ 
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Remark 10. In the above theorem we have P n = , because qj = if I ^ 0. Then 

[Pj,P n ] = implies [Pj, A {lj ... |n} ] = 0. 

Hierarchy of Integrable Potentials of Type D n (n > 3) 

Toda-P>i^ — ► Toda-.D n 

Ellip-Dn Trig-D n -> Rat-£> n 

7 Classical limits 

For functions f(£,x) and <?(£, x) of (£, x) = (£i, . . . , £ n , xi, . . . , x n ), we define their Poisson 
bracket by 



fc=i 

Theorem 13. Put 



df_dg_ _ dg df 

d£k dx k d£ k dx k 



1 n 

P{M = --Y,e k +R{x). 



2 

k=l 

Then for the integrable potential function R{x) given in this note, the functions P k {^,x) and 
A{i,...,n}(£,i x ) of (£,,x) defined by replacing d v by £ u (y = 1, . . . , n) in the definitions of P k and 
A{i,...,n} in Sections [3J [7] an d satisfy 

{P i (Z,x),P j (Z,x)} = {P(Z,x),P k (Z,x)}=0 for \<i<j<n and 1 < k < n. 

Hence P(£, x) are Hamiltonians of completely integrable dynamical systems. 
Moreover if the potential function R(x) is of type D n , then 

{A {1) ... )n} (£,z),P fc = {A {1 ^ n} (Z,x),P(Z,x)} =0 for l<k<n. 

Proof. If R(x) is a potential function of (Ellip- J 4 n _i), (Ellip-i? n ) or (Ellip-D n ), the claim is 
proved in |29[ [32] . Since the claim keeps valid under suitable holomorphic continuations with 
respect to the parameters which are given in the former sections, we have the theorem. ■ 

Remark 11. Since our operators P k are expressed by operators P£ = X^Pfc i( x )q ki (d) such that 

i 

the polynomials q ki (d) satisfy \p k i{x),q ki {d)} = 0, there is no ambiguity in the definition of 
the classical limits by replacing d u by £ u . In another word, if we have given the above integrals 
Pj(x,£) of the classical limit, we have a natural unique quantization of them. 

8 Analogue for one variable 

Putting n = 1 for the Schrodinger operator P of type A n in Section [3] or of type B n in Section [5J 
we examine the ordinary differential equation Pu = Cu with C 6 C (cf. [411 § 10.6]). We will 
write the operators Q = P — C . 

(Ellip-fii) The Heun equation (cf. [32l § 8], gH pp. 576]) 

Id 2 3 
"2^2 +£C>(t + u; J )-C. 

i=o 
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(Ellip-Ai) The Lame equation 

1 d 2 t , , ^ 

(Trig-i?Ci-reg) The Gauss hyper geometric equation 

2 <it 2 sinh 2 At sinh 2 2At 
(Trig-Ai) The Legendre equation 

Id 2 Co 
2cft 2 sinh 2 Ai 

(Trig-Bi) with C = Ci = C 3 = 0. The (Modified) Mathieu equation 

1 d 2 

--— ? + C 2 cosh2At-C. 

2 ar z 

(Rat-Si-2) Equation of the paraboloid of revolution 

1 d 2 Co _ 2 _ 
"2^ + ^ + Clt " C - 

This is the Weber equation if Co = 0. With s = t 2 and using the unknown function t^u, the 
above equation is reduced to the Whittaker equation: 

| C'q C[ 

2 ds 2 s 2 s 

(Rat-ylo-bry2) with C 2 = C 3 = 0: 

-a 

If C\ 7^ 0, this is transformed into the Weber equation under the coordinate s = t + C2/(2C\). 
If C\ = 0, this is the Stokes equation which is reduced to the Bessel equation. In particular, the 
Airy equation corresponds to C = C\ = 0. 
(Toda-.BCi) 

which is transformed into (Rat-2?i-2) by putting s = e~ l . In particular 
(Tbda-Ai) 

is reduced to the Bessel equation. 
(Rat-Ai) the Bessel equation 
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In fact, the equation —u"/2 + Cou/t 2 = Cu is equivalent to 



since 



_i d 2 i d 2 Id 1 

1 2 °d^ oU ~ ~d£ 2+ Ut~ W 

Hence if C ^ 0, the function v = t^' 2 u satisfies the following Bessel equation with s = \/—2Ct: 

d 2 v | ldv / C + l/8 ^ 
ds 2 s ds \ Cs 2 

Hierarchy of ordinary differential equations 

Heun => Lame — > Legendre 

\ fo:2 \ 

2'5 5 - 3 2'1 

Mathieu <= Trig-£>i Gauss Bessel (Stokes) — > Airy 



|5:4 \ \ lb:2 h: 



1 



Rat-Ao-bry Rat-Bi U Whittaker ^1 Weber 



9 A classification 

We present a conjecture which characterizes the systems listed in this note. 

Let P be the Schrddinger operator with the expression (jl.ip and consider the condition 



there exist P±, . . . ,P n such that < 



p g c[Pi, . 



i -F71] i 



[1 < i < j < n) 



(9.1) 



y e 2 ■■■e 2 

^ s ji 

l<ji<-<3k<n 



(1 < k < n). 



Note that all the completely integrable systems given in Sections [31 H] or [5] satisfy this condition. 
Conjecture. Suppose P satisfies (19. ip . Under a suitable affine transformation of the coordinate 



x G C n which keeps the algebra 



ydiyd£,...,yd 2 



2/i 



invariant, P is transformed into an 



L k=i k=l k=l 

integrable Schrddinger operator studied in Sections{3l\4\ or\^ (namely ufj, and Wf. in (jl.3p are 
suitable analytic continuations of the corresponding functions of the invariant elliptic systems) 
or in general a direct sum of such operators and/or trivial operators 



d 2 



with arbitrary functions v{x) of one variable. 

Here the direct sum of the two operators Pj(x,d x ) = y a a (x)d x of x G C Uj for 

j = 1, 2 means the operator Pi(x, d x ) + P2{y, d y ) of (x, y) G c ni+n2 . 

We review known conditions assuring this conjecture and give another condition (cf. Theo- 
rem [TO] and Remark I17|) . We also review related results on the classification of completely 
integrable quantum systems associated with classical root systems. 
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Remark 12. The condition 

there exists P 2 such that [P,P 2 ] = and a(P 2 ) = ^ (9.2) 

l<i<j<n 

may be sufficient to assure the claim of the conjecture. 
Remark 13 (Type A2). If n = 2 and if there exists P3 satisfying 

= 66 + 66 + 66 and [p,p 3 ] = [ai + a 2 + a 3 ,p] = + d 2 + a 3 ,p 3 ] = o, 

then Conjecture is true. 

In fact this case is reduced to solving the equation 



u(x) u'{x) 1 
v(y) v'(y) 1 
w(z) w'{z) 1 



for x+y+z=0 (9.3) 



for three unknown functions u(t), v(t) and w(t), which is solved by [311]. Here u(t) = u ei _ e2 (t), 
v(t) = u e2 - e:j (t) and 10 (i) = n ei _ e3 (-t). 

9.1 Pairwise interactions and meromorphy 

Theorem 14 ( [40] ) . The potential function R(x) of P satisfying (|9.ip zs 0/ i/ie /orm 

P(x) = Yl (9.4) 

a£S(B„)+ 

w'i/i meromorphic functions u a (t) of one variable. 
Remark 14. i) The condition (19.21) assures 



„2 +7^2 J +^(si,x 2 ) 



■[XjXfc 

a£Y,(B n )+ l<i<j<k<n 

with Cijk £ C and thus the above theorem is proved in the invariant case (cf. Section [H2]) by 
or in the case of Type P 2 by [23] or in the case of Type A n —\. This theorem is proved in [30 
by using [P,P 2 ] = [P,P 3 }=0. 

ii) Suppose n = 2 and the operators 

P - A(&_ I 2 

2 \dx1 dx'2 

nn 9 m m , \ ^ 

^ rir l rir m ~ l ^ rl^ Pit 3 

satisfy [P,T] = and <r m (T) ^ C[er(P)]. Here c,- G C. Then (SB Theorem 8.1] shows that there 
exist functions u u ^{t) of one variable such that 

L m„ — 1 

R(xi,x 2 ) = ^ (&„aii + a^^u^a^i - 6^x 2 ) 
by putting 



% - ^ e ^ m ~ v = n m - ^r- • 

V W i=0 v=l 

Here (a„, 6^) 6 C 2 \ {(0, 0)} and a u b^ ^ a^b u if \l ^ v. 
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Definition 6. By the expression (|9.4p . put 

S = {a 6 E(P„)+ ; u' a ^ 0} 
and let W(S) be the Weyl group generated by {w a ; a G 5} and moreover put 5 = W(/S)<S. 

Theorem 15 ([23] for Type B2, [10] in general). If the root system S has no irreducible compo- 
nent of rank one, then (|9.2p assures that any function u a (t) extends to a meromorphic function 
on C. 

Remark 15 ([32] (6.4)-(6.5)], 03 § 3]). The condition flOD is equivalent to 

S« = S 3 * (1 < i < j < n) (9.5) 

with 

S % 3 = \dfvi{xi) + d l { u t( X i + x v) + U iv( X i ~ x »)) i. u tj( x i + x i) - u ij( x i ~ x j)) 

V vei(ij) / 

+ 3 \diVi(xi) + ^ di(uf u (xi + x v ) + u~ v {xi - x v ))\ (diufjfa + xj) - diU~j(xi - Xj)) 



+ 2lt; 4 (x i )+ J] {ul{ Xi + x v ) + 

V wG/(ij) 



- ( d i u tj(. x i + - d i u ij( x i ~ x j)) 



+ Yl ( d i u tu( x i + x u)-dfu-(x i -x u ))(uj u (x j + x u )-uj u (x j -x u )). 
vei(i,j) 

Here j) = {1,2, ... ,n}\{i,j}. 

Lemma 5. Suppose P satisfies (|9.2p and (|9.4p . Lei So 6e a subset of S such that 



S C Y Re i and S\S C Y Re '' 

i=l i=m+l 

with a suitable m. Then the Schrodinger operator 

1 m 

i=i ae5 ns 

m 

on M m admits a differential operator P' 2 on M m satisfying \P' , P 2 ] = and o~{P!^) = Yl £i£j> 

l<i<j<n 

that is, the condition (19. 2|) uwi/i replacing P by P' . 

Proof. This lemma clearly follows from the equivalent condition (I9.5P given in Remark [15] ■ 



9.2 Invariant case 

Theorem 16 ( [24] l25l [29] 132]). Assume that P in (II. ip is invariant under the Weyl group 
W = W{A n _i), W{B n ) or W{D n ) with n > 3, or = W(B 2 ). // we have (T2J) twift 



Pi = ft + d 2 + ■ ■ ■ + d n if W = W(A n _i), 
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<r(Pk) 



£ if W = W{A n _ 1 ) and 1 < k < n, 

l<ii<J2<---<jfc<n 

E ^ W = W(B n ) and 1 < k < n, 

l<7l<J2<---<jfc<n 

E */ W = W(£» n ) and 1 < A; < n, 

U<jl<j2<---<jfc<n 



<7(Pn)=ei6---en </ W = W(A0, 
Conjecture is true. 
Remark 16. The condition 

[P, Pi] = [P, P 3 ] = if W = W(A n ^), 

[P, P 2 }=0 if W = W{B n ) or W = W{D n ) 

together with (|9.4|) is sufficient for the proof of this theorem. 

9.3 Enough singularities 

Put H = {a£ Yi(B n ) + ; u a (t) is not entire}. 

Theorem 17. i) (|23j) Suppose n = 2 and let S be of type B>2- If #H > 2, i/ien Conjecture is 
true. 

ii) ([10]) If S is of type A n _\ or of type B n and moreover the reflections w a for a £ S generate 
W{A n _i) or W(B n ), respectively, then Conjecture is true. 

This theorem follows from the following key Lemma. 

Lemma 6 (|23 t 137} I40j). Suppose (|9.1|) and moreover that there exist a and (3 in S such that 
a /3, {a, (3) ^ and u a (t) has a singularity at t = to. Then u a (t — to) is an even function 
with a pole of order two at the origin and 

u Waif3) (t - 2to|^y) = u p {t) if w a {f3) E £(B n )+ 

u_ WaW (-t + 2t o-|^4) = u p (t) if ~ W M G S(^n) + . (9.6) 

Corollary 1. Suppose the assumption in Lemma\^ 

i) If u a {t) has another singularity at t\ ^ to, then 

uJ t + 2(h -to)p^\ = Ul {t) for 7 eS. (9.7) 
V {a, a) J 

ii) Assume that u a has poles at 0, to an d t\ such that to and t\ are linearly independent 
over K. Then up{t) is a doubly periodic function and therefore up(t) has poles and hence u a (t) 
is also a doubly periodic function. We may moreover assume that up has a pole at by a parallel 
transformation of the variable x. 

Case I: Suppose a = ei — ej, (3 = ej — with \ <i<j<k<n. 

u ei - ej (t) =U ej -e k (t) =U ei -e k (t) = Cp{t;2u\,2u)2) + C 

with suitable C, C E C, which corresponds to (Ellip-^)- 

Case II: Suppose a = — ej and (3 = ej with 1 < i < j < n. 

Then ( Uei - ej (t),u ei+ej (x),u ei (t),u ej (t)) is (Ellip-J^), (Ellip-^-S) or (Ellip d -S 2 ). 
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iii) If S is of type A n _i or B n or D n and one of u a (t) is a doubly periodic function with 
poles, then P transforms into (Elhp-A n _x) or (Ellip-£? n ) or (Ellip-.D n ) under a suitable parallel 
transformation on C n . 

Proof of Corollary Q3 i) is a direct consequence of Lemma [6j iii) follows from ii). We have 
only to show ii). 

Case I: It follows from (|9.6p that u a (t) = up{t) = u ei _ ek (t) and they are even functions. Let 

^2u u 2u 2 = {2miuji + 2m 2 w 2 ; mi,m 2 G Z} 

be the set of poles of u a . Then (|9.7p implies up(t + 2lo±) = up{t + 2uj2) = up{t). Since 2u\ and 
2lo 2 are periods of p(t) and there exists only one double pole in the fundamental domain defined 
by these periods, we have the claim. 

Case II: It follows from (j9.6j) that u ei _ ej (i) = u ei+ej (t) and u ei (t) = u ek (t) and they are even 
functions. Let T 2 u u 2ui 2 be the poles of u ei - ej (t). Then (|9.7|) means u ei (t + 2uii) = u ei {t + 2uj 2 ) = 
u ei (t). Considering the poles of u ei - ek (t) with (|9.7p . we have four possibilities of poles of u ei : 

(Case II-O): r 2wi ,2a, 2 , 
(Case II-l): T 2u)1 , 2u)2 LK^i + ^2wi,2w 2 ), 
(Case II-2): T 2u)1 , 2u)2 U( w 2 + ^2w u 2w 2 ), 
(Case II-3): T 2u)1 , 2u)2 LK^i + ^1,2^2) U( w 2 + I^wi,^)- 
Here we note that (Case II- 1) changes into (Case II-2) if we exchange uj\ and u> 2 . Then we have 

(Case II-O): u £i _ ej (t + 4wi) = u ei _ e ,.(i + 4w 2 ) = u(t), 
(Case II-2): u ei _ ej (t + 4wi) = u ei - ej (t + 2w 2 ) = tt(t), 
(Case II-3): u e ._ e .(t + 2u>i) = u ei _ e .(t + 2u 2 ) = u(t). 

Thus (Case II-O), (Case II-2) and (Case II-3) are reduced to (Ellip d -£ 2 ), (Ellip-^-S) and 
(Ellip-£? 2 ), respectively. ■ 

Let TL be a finite set of mutually non-parallel vectors in W 1 and suppose 

p = -\ E 5 I + = £ ^fel + 

Here C a are nonzero complex numbers and R(x) is real analytic at the origin. We assume that 
TL is irreducible, namely, 

R" = ^ Ma, 

0/V7f'^^3aeH' and 3[3eH\H' with (a, 0)^0. 

Definition 7. The potential function _R(x) of a Schrodinger operator is reducible if i?(x) and 
M n is decomposed as R(x) = R\(x) + R 2 (x) and M n = V\ © V" 2 such that 

0Cy 1 C]R n , V r 2 = V r i L , d V2 Ri{x) = d Vl R 2 (x) = for Vu 2 G F 2 and Vi>i G Vi. 

If i?(x) is not reducible, R(x) is called to be irreducible. 

Theorem 18 ([37]). Suppose n > 2 and there exists a differential operator Q with [P,Q] = 
whose principal symbol does not depend on x and is not a polynomial of ^ £?. Put W = 
{w a ;a G TL}. If 

2C a / k(k + 1) for keZ and aeTL, (9.8) 
then W is a finite reflection group and a{Q) is W -invariant. 
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9.4 Periodic potentials 

The following theorem is a generalization of the result in [30] . 

Theorem 19. Assume R(x) is of the form (|9.4p with meromorphic functions u a (t) on C and 
/ 2ir\/—la \ 

R[x+ , = R{x) for a G E(B n ) (9.9) 



(a, a) 

and moreover assume that 

the root system S does not contain an irreducible component of 

type B2 or even if S contains an irreducible component 

S2 = {±ej ± ej, ±ej, ie^} of type B2, the origin s = is not an isolated 

essential singularity of u a (log s) for a 6 ^f] X(D n ) + . 

Then Conjecture is true. 

Remark 17. i) The integrable systems classified in this note which satisfy the assumption of 
Theorem [19] under a suitable coordinate system are (Ellip-*) and (Trig-*) and (Toda-*), which 
are the systems given in this note whose potential functions are not rational. 

ii) The assumption (|9.9p implies that u a (logs) is a meromorphic function on C \ {0} for any 
a G T,(B n ) + . It means that the corresponding Schrodinger operator is naturally defined on the 
Cartan subgroup of Sp(n, C) with a meromorphic potential function. 

Lemma 7. Assume n = 2, (|9.2p . (j9.9j) . S is of type B2 and moreover n Q (logs) are holomorphic 
for a £ X(£?2) + and < |s| <C 1. // the origin is at most a pole of upilogs) for (3 E S(Z?2) + ; 
the origin is also at most a pole of u a (logs) for a G X(£?2) + - 

Proof. Use the notation as in (14.21). Put 



u 



(log s) =Uq+J2 vU-s v , u+(log s) = U++J2 

v=r v=m 

00 00 

v(logs) = V + uV " sU > w(logs)=W + vW u s u . 

v=—oo v=—oo 

with U~ , , V u , W u G C, rm / and (U~, C/+) / 0. Then as is shown in [30] the condition 
for the existence of T(x, y) in (|4.3p is equivalent to 

pq(2p -q)(p- q)(V 2p - q U+_ p + V q U~_ q + W q - 2p U+ - W q U~) =0 for p, q G Z. (9.10) 

Hence if p < r and p < m, 

P(p ~k)(p + k)k(V p+k U+ k + V p - k U k ) = for k£Z. 

Case U~ 7^ 0: Put k = r. Suppose q is negative with a sufficiently large absolute value. 

00 

Then V q = (—U^ r /U~)V q+ 2r, which implies V q = since vV u s' / converges for < \s\ <S 1. 

v=~oo 

Suppose q is negative with a sufficiently large absolute value compared to p. Then by the 
relation W q -2 P U£ — W q U~ = we similarly conclude W q = 0. 

Case ^ 0: Putting k = —m, we have the same conclusion as above in the same way. ■ 
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Proof of Theorem 1191 Lemma [5] assures that we may assume S is an irreducible root system. 
We may moreover assume that the rank of S is greater than one. 

Suppose that there exists 7 G S such that the origin is neither a removable singularity nor 
an isolated singularity of u 7 (log s). Then Uy(t) is a doubly periodic function with poles. Owing 
to Corollary [H ]P u a ((a,x)) is reduced to the potential function of (Ellip-A n _i) or (Ellip-P n ) 

or (Ellip-A,). 

Thus we may assume that the origin is a removable singularity or an isolated singularity of 
■Uq, (logs) for any a G S. 

Let a, (3 G Sf]T,(D n ) with a ^ (3 and (a, (3) 7^ 0. Put 7 = w a f3 or 7 = —w a (3 so that 
7 G £(-D n ) + . Then [30] shows that u(t) = u a (t), v(t) = up(t) and w(t) = u 7 (— i) satisfy (|9.3j) . 
Then Remark 1131 savs that the origin is at most a pole of u(logs), f;(logs) and w(logs). 

Let a G Sf]T,(D n ) and (3 G S \ T,(D n ) with (a,/3) / 0. Let W be the reflection group 
generated by w a and wp and put 5° = VF{a, (3} f) T,(B n ). Then [30] also shows that 

defines an integrable potential function of type B2. Hence Lemma [7] assures that the origin is 
at most a pole of u a (\og s) for a G 5°. 

Since S is irreducible, the origin is at most a pole of ii a (logs) for a G S. Then Theorem [T9l 
follows from [30]. ■ 

9.5 Uniqueness 

We give some remarks on the operator which commutes with the Schrddinger operator P. 
Remark 18 ([32, Lemma 3.1 h)]). If differential operators Q and Q' satisfy [Q, Q'\ = 0, cr(Q') = 

n 

£,f and ord(Q) < N — 2, then Q has a constant principal symbol, that is, a(Q) does not 

3=1 

depend on x. 

Hence if there exist differential operators Qi, . . . ,Q n with constant principal symbols such 
that a(Qi), . . . , a(Q n ) are algebraically independent and moreover they satisfy [Qi, Qj] = for 
1 < i < j < n, then any operator Q satisfying [Q,Qj] = for j = 1,... ,n has a constant 
principal symbol. In particular, if a differential operator Q satisfies [Q,-Pfc] = for in (jl.2p 
and p.4p with k = 1, . . . ,n, then a(Q) does not depend on x. 

Remark 19. Assume that a differential operator Q commutes with a Schrddinger operator P 
and moreover assume that there exist linearly independent vectors Cj G C n for j = 1, . . . , n such 
that the operator is invariant under the parallel transformations x >— > x + Cj for j = 1, ... ,71. 
Then a{Q) does not depend on x (cf. [32j Lemma 3.1 i)]). 

Furthermore assume that P is of type (Ellip-i 7 ) or (Trig-i 7 ) or (Rat--F) with F = A n _\ or B n 
or D n . If the condition (|9.8p holds or Q is VF(i ? )-invariant, it follows from Theorem 1181 or [291 
Proposition 3.6] that Q is in the ring C[Pi, . . . , P n ] generated by the VK(i ? )-invariant commuting 
differential operators. If the condition (8.11) is not valid, cr(Q) is not necessarily W(i ? )-invariant 

(cf. [2E51CS]). 

Remark 20 ([32, Theorem 3.2]). Let P be the Schrddinger operator in Theorem[l6j Under the 
notation in Theorem 1161 suppose P^ are VK-invariant for 1 < k < n. Then the ring C[Pi, . . . , P n ] 
is uniquely determined by P and Q, where Q = P3 if W = W(A n -i) and Q = P2 if W = W(B n ) 
or W(D n ). 
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Remark 21. If P c = —(1/2) ^ d 2 + cR(x) is a Schrodinger operator with a coupling constant 

7=1 

c G C such that P c admits a non-trivial commuting differential operator Q c of order four for any 
c G C, then the operator P c may be a system stated in Conjecture under a suitable coordinate 
system. 

The following example satisfies neither this condition nor the condition (|9,8p . It does not 
admit commuting differential operators (|1.2p satisfying (|1.4p if m ^ 0, —1. 

Example 2. It is shown in [BJ [35] that the Schrodinger operator 

-i n n2 / i -i \ n ~l , 1 

l y\ a m(m + lj m + 1 

is completely integrable for any m and algebraically integrable if m is an integer. 

The following example shows that the Schrodinger operator P does not necessarily determine 
the commuting system C[P\, . . . , P n ]. 

Example 3. Let a, fl, 7 and A be complex numbers. Put (Aq, A±, Co, C\) = (a, 7/2 — A/2, (3, A) 
for (Rat-i?2-S) in Theorem [3] (cf. [32 \ Remark 3.7]). Then the Schrodinger operator 

Fa ^~~2\bW + dy^J +{X +V \(x 2 -y 2 ) 2 + x^ 2+1 
commutes with 

( 92 4aX V n/ n Y J P x l\ 92 J P x 2\ 92 



Q «to* = {m + o^w " 2(7 " x)xy ) " 2 (? + v ) w 2 (^-- Vri 

^2 + Xx ){^ + Xy ) + — (x 2 _ ^2)2 — + 8A (^ - x ) x y 



dy 



2 



for any A G C. Note that [Q a ,/?,7,Ai Qa,/3,7,A'] / if A 7^ A' and these operators are W(B2)- 
invariant. The half of the coefficient of the term A of <3 Qi/ 3 i7i a considered as a polynomial 
function of A is 

( d d \ 2 ( xy xy \ I ' y 2 x 2 \ 22 

Sa ^ = -{ y d^- x dy-J + 2a \V^W " ¥TW ) +w \~ 2 + V') +Alxy ■ 

In particular, P = — (l/2)(<9 2 + d 2 ) + 7(x 2 + y 2 ) commutes with d x d y — 2~/xy and xd y — yd x . 

- n 

Note that if R(x) is a polynomial function on C n , the condition —(1/2) ]P <9 2 + /?(#), Q 

for a differential operator Q implies that the coefficients of Q are polynomial functions (cf. 
Lemma 3.4]). 

9.6 Regular singularities 

Definition 8 ([IB]). Put tf fc = t k d/dt k and Y" fc = {t = (*i,...,t„) G C n ; i ft = 0}. Then 
a differential operator Q of the variable t is said to have regular singularities along the set of 
walls {Y u ...,Y n } if 



Q = g(0i,...,0„) + 5NfcQfcM). 



fc=i 
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Here q is a polynomial of n variables and Qt are differential operators with the form 

and a a (t) are analytic at t = 0. In this case we define 

cr*(Q) = ...,£„) 
and cr*(Q) is called the indicial polynomial of Q. 

Theorem 20. Ze£ R(t) be a holomorphic function defined on a neighborhood of the origin ofC n . 
Let Qi and Q2 be differential operators oft which have regular singularities along the set of walls 
{Yi, . . . , Y n }. Suppose a*(Qi) = c*(Q2) and [Q\, P] = [Q2, P] =0 with the Schrddinger operator 

Then Q\ = Q 2 . 

Proof. Put tj = e~( Xj ~ Xj+1 ^ for j = 1, . . . ,n — 1 and t n = e~ Xn . Then dj = — i)j for j = 
1, ... ,7i — 1 and d n = — i? n . Under the coordinate system x = (x%, . . . , x n ) Remark 1 191 savs that 
Qi — Q2 has a constant principal symbol, which implies Q\ = Q2 because ct*(Qi — Q2) = 0. ■ 

A more general result than this theorem is given in [31j . The following corollary is a direct 
consequence of this theorem. 

Corollary 2. Put tj = e~ x ( Xj ~ Xj+1 ' for j = 1, . . . ,n — 1 and t n = e~ Xxn . Suppose P is 
the Schrddinger operator of type (Trig-j4 n _i), (T>ig-j4 n _i-bry-reg), (Trig-i?C n -reg), (Trig-.D n ), 
(Toda-A n _i), (Toda-£C n ) or (Toda-D n ). 

i) P and P/. for k = 1, . . . , n have regular singularities along the set of walls {Y%, . . . , Y n }. 

ii) Let Q be a differential operator which has regular singularities along the set of walls 
{Yi, . . . , Yn} and satisfies [Q, P] = 0. If a*(Q) = <r*(Q) for an operator Q G C[Pi, . . . , P n ], then 
Q = Q. 

Remark 22. i) This corollary assures that certain radial parts of invariant differential operators 
on a symmetric space correspond to our completely integrable systems with regular singularities 
and the map cr* corresponds to the Harish-Chandra isomorphism (cf. [31]). 

ii) The system (Trig-i?C n -reg) is Heckman-Opdam's hypergeometric system [11] of type BC n . 
Since (Trig-i?C n -reg) is a generalization of Gauss hypergeometric system related to the root 
system S(B n ), the systems in the following diagram are considered to be generalizations of Gauss 
hypergeometric system and its limits (cf. Section [8]). They form a class whose eigenfunctions 
should be easier to be analyzed than those of other systems in this note. 

Hierarchy starting from (Trig-BC^-reg) 



Toda-A n _i 

/ 

Rat-D ?t Trig-^ n _i -> Rat-A n _i 

ft3:l 1>3:1 1>3:1 

Rat-5 n -2 Trig-^ n _i-bry-reg -> Rat-,4 n _i-bry2 

T / 

Trig-£C n -reg -> Toda-D n -bry -» Toda-5C n 

JJ-3:1 JJ-3:1 JJ-3:1 

Trig-D n — > Toda-D n — > Toda-A n _i 

1 \ 
Trig-A n _! Rat-D n 
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9.7 Other forms 

If a Schrodinger operator P is in the commutative algebra D = C[Pi, . . . , P n ], then the differential 
operator P := Tp(x)~ 1 P o ^(x) with a function i s hi the commutative algebra D = 

C[ip(x)~ 1 P± o ip(x), . . . ,'0(x)" 1 P n o of differential operators. Then 



^ = "a E + E^) + R ^ ( 9 - u ) 

3=1 •? 3=1 3 

dip(x 



dxj 



dj(x) for j = l,...,n. (9.12) 



Conversely, if a function ip(x) satisfies (|9.12p for a differential operator P of the form (|9.1ip . 
then P = iJ){x)P o ^(x) -1 is of the form (jl.ip . which we have studied in this note. 
If i/j(x) is a function satisfying 

2ip(x) ^ dx\ 
j — i ^ 

then 



P = V0*0 



1 f 4 E 5 * + R w) ^) = -5 E 5 f - ^r 1 E 
\ j i / i=i j=i ^ 



Note that 



dxj dxj ' ^— ' cte 2 ' cAr 2 ' V cte 7 - 

j j j=i j j=i j jf=i v j 



Putting 



(x) = m logsinhA(xj — Xj) 

l<i<j<n 



we have 



Am cothA(xfc 



dxi 

l<i<n, ij^k 

n »2 



E^+E(l?) 2 -^ E 



j=l 3 fe=l * \<i<j<n 



+ 2A 2 m 2 ^ coth 2 X(xi - xj) + A 



2 2 ~ !)( n - 2 ) 



3 

l<£<j<n 

= 2A 2 m(ra-l) sinh" 2 A(ccj - Xj) + A 2 m 2 — - 

3 

l<i<j<n 

since 

coth a ■ coth /3 + coth (3 ■ coth 7 + coth 7 • coth a = — 1 if a + /? + 7 = 0. 
Hence 



1 ™ 

P = — -^^(9 2 — m A coth A(xj — Xj)(di — djj, 



2 

3=1 l<*<3f^ n 
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ij){x)= A m sinh m A(xi-Xj), 

l<i<j<n 

, i / x 1 \— > «9 mini — 1)A 2 m 2 n(n 2 — 1)A 2 



and P is transformed into the Schrddinger operator of type (Trig-^4 n _i). 
Now we put 

<j)(x) = itiq (log sinh A(x» — Xj) + log sinh A(x« + Xj)) + m\ log sinh Ax^ 

l<i<j<n l<k<n 



+ ni2 ^] log sinh 2Axfc 



Kfe<n 



and we have 
d(j)(x) 



- Xmo (coth A(xfc + Xj) + coth A(x& — Xj)) + Ami coth Ax^ + 2Am,2 coth 2Ax/j, 

coth Ax& coth 2Axfc = 1 + — sinh" 2 Xxk, 

(2 coth A(x& + Xj) coth A(x& — Xj) + 2 coth A(x& + Xj) coth A(x& — xj) 

{i,j,k}=I 

+ coth A(xfc + Xj) coth A(xfc + Xj) + coth A(xfc — Xj) coth A(x^ — Xj)) 

= (coth A(xfc + Xi) coth A(xfc + Xj) + coth A(xj — Xj) coth A(xj + Xfc) 

+ coth A(xj — Xj) coth A(xj + Xk)) + (coth A(x& — Xj) coth A(x& — Xj)) = 8 

{i,3,k}=I 

for 7c{l,...,n} with # J = 3, 

cothA(x fc + Xj) + cothA(x fc - Xj) = smh2Ax fc 

smn A(Xfc + Xi) sinh A(x^ — Xj) 

cosh 2Axfc — cosh 2Axj 2 cosh 2 Ax^ — 2 cosh 2 Axj 



sinh A(xfc + Xj) sinh A(x/c — Xj) sinh A(x& + Xj) sinh A(x& — Xj) 
E(~1^~0 = 2A2?71 o X] (coth 2 A(xj - Xj) + coth 2 A(xi + Xj)) 

fc=l ^ fe ' l<i<j<n 

n n n 

+ \ 2 m\ coth 2 Ax?; + 4A 2 77i2 coth 2 2Axfe + c l\ 2 m\vti2 sinh" 2 Axfe 

fc=i fc=i fc=i 

4X 2 m 2 ] n(n — l)(n — 2) , n „ 
H u v — + 2A^m (mi + 2m 2 )n(n - 1) + 4\ 2 m 1 m 2 n, 

E ^7T^ +Ef^V= 2A 2 mo(m - 1) £ (sinh" 2 A(x,-Xj) + sinh" 2 A(x J +x J )) 

3=1 6X 3 fe=l V ° Xk J l<i<j<n 

n n 

+ A 2 mi(mi + 2?tt-2 — 1) sinh -2 Ax/c + 4A 2 m2(m,2 — 1) sinh -2 2Axfc 

fc=i fc=i 

+ A 2 ( ( — mo(2n — 1) + 2mi + 4m2 )mo(n — 1) + (mi + 2m2)\n. 
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Hence 
P 



_. n n / 

o ^ &j - ^ J A I ^ mo(cothA(xi-Xfc)+cothA(xi + x fc )) 

j=l fc=l \l<i<j<n 



+ mi coth Axfc + 2m2 coth 2Axfc J <9fc, 

n n 

i>(x) = (sinh m ° X( Xi - xj) sinh m ° A(x* + xj)) JJ sinh mi Ax fc JJ sinh™ 2 2Ax fc (9.14) 

l<j<j<n fc=l k=l 

and P is transformed into the Schrodinger operator of type (Trig-I?C n -reg): 

1 n /A 2 A 2 \ 

v - mi(mi + 2rri2 — 1)A 2 \ ^ 2m2(m2 — 1)A 2 
2sinh 2 Ax fc sinh 2 2Ax fc 

fe=l k=l 



+ A 2 (^°(2n - l)(n - 1) + m (mi + 2m 2 )(n - 1) + (mi + 2m2 ^ 



n. 



Remark 23. As is shown in [151 Theorem 5.24 in Ch. II], the operator (I9.13P or (19. 14ft gives the 
radial part of the differential equation satisfied by the zonal spherical function of a Riemannian 
symmetric space G/K of the non-compact type which corresponds to the Laplace-Beltrami 
operator on G/K. Here G is a real connected semisimple Lie group with a finite center, K is 
a maximal compact subgroup of G and the numbers 2m, 2mo, 2mi and 2m2 correspond to the 
multiplicities of the roots of the restricted root system for G. 

Similarly the following operator P is used to characterize the i^-fixed Whittaker vector v on 
G = GL(n,R) 

1 ™ JL /r> 4- 1 i \ S-T , * E (j/2-n+l/4)xi 



^ = - ^ E 5 f + E ( ^ - 1 ) + c E ^ > = ei 



2^ J ^ 48 

i=i i=i 

Namely v is a simultaneous eigenfunction of the invariant differential operators on G/K and 
satisfies v(nx) = x( n ) v ( x ) with n £ N and x G G/K. Here G = KAN is an Iwasawa de- 
composition of G and x is a nonsingular character of the nilpotent Lie group N. Then v\a 
is a simultaneous eigenfunction of the commuting algebra of differential operators determined 
by P. 
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